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ABSTRACT 

We used the Submillimeter Array (SMA) to image 860 /im continuum and CO (3-2) line emission in 
the ultraluminous merging galaxy Arp 220, achieving a resolution of 0'.'23 (80 pc) for the continuum 
and 0'.'33 (120 pc) for the line. The CO emission peaks around the two merger nuclei with a velocity 
signature of gas rotation around each nucleus, and is also detected in a kpc-size disk encompassing the 
binary nucleus. The dust continuum, in contrast, is mostly from the two nuclei. The beam-averaged 
brightness temperature of both line and continuum emission exceeds 50 K at and around the nuclei, 
revealing the presence of warm molecular gas and dust. The dust emission morphologically agrees 
with the distribution of radio supernova features in the east nucleus, as expected when a starburst 
heats the nucleus. In the brighter west nucleus, however, the submillimeter dust emission is more 
compact than the supernova distribution. The 860 /im core, after deconvolution, has a size of 50-80 
pc, consistent with recent 1.3 mm observations, and a peak brightness temperature of (0.9-1.6) x 10 2 
K. Its bolometric luminosity is at least 2 x 10 11 L & and could be ~10 12 i© depending on source 
structure and 860 /im opacity, which we estimate to be of the order of T86o ~ 1 (Le., Nn 2 ~ 10 25 
cm -2 ). The starbursting west nuclear disk must have in its center a dust enshrouded AGN or a very 
young starburst equivalent to hundreds of super star clusters. Further spatial mapping of bolometric 
luminosity through submillimeter imaging is a promising way to identify the heavily obscured heating 
sources in Arp 220 and other luminous infrared galaxies. 

Subject headings: galaxies: ISM — galaxies: starburst — quasars: general — galaxies: evolution — 
galaxies: individual (Arp 220) 



1. INTRODUCTION 

Arp 220 is the nearest ultraluminous infrared galaxy 
(ULIRG) and is an advanced merger (I/s_iooo nm — 
1O 12 2 L ; D = 75 Mpc; 1"= 361 pc; see Table P). Its 
vast luminosity is radiated almost entirely between in- 
frared and millimeter wavelengths as thermal dust emis- 
sion. This is because the source, an extreme starburst, 
a quasar-class AGN(s), or both, is deeply embedded 
in dust within the central few kpc (|Soifer et alj 119841 : 
I.Tov et al.lll98l . This merger has been a prime target to 
study the mechanism to generate the high luminosity of 
ULIRGs and to study galaxy evolution through majo r 
mergers (|Sanders fc Mirabell Il996t iGenzel et al. l fl998h . 
Detailed studies of Arp 220 as well as other local ULIRGs 
may also help us to understand submilli meter-detected 
lumin ous galaxies at high redshift (e.g., ICoppin et al.l 
120061 ). beca use the latter may be scaled-up versions of 
the former (jTacconi et al.ll2006] ). 
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The Arp 220 system has two nuclei, presumably 
from the merger progenitors , with a projected separa- 
tion of about 350 pc (|Norrislll988T: iGraham et alj|1990t 
iClements et al.ll2002[ ). Millimeter CO observations reveal 
a large amount of molecular gas (several 1 9 in the 
central few kpc with pea k s at the nuclei (|Scoville et al 



fl99l ]lScoville et al.lll997l fdSYB97ft: lD ownes fc Solomon! 
mm Sakamo toet al.l 119991 fdSS99): D ownes fc Eckart 



120071 ( dDE07D ). Sub-arcsecond resolution CO and H I 
observations show that a small gas disk (r ~ 50 pc) ro- 
tates around each nucleus and that the two nuclear disks, 
counterrotating with respect to each other, are encom- 
p assed by a larger outer disk (|SS99l : iMundell et"ai1l2001l 
( dMFPOlf ), in which one of the nuclear disks was mod- 
eled as a part of the outer dis k.) Th e two nuclei dominate 
the dus t emission at 1.3 m m (lSS99f) and in mid-IR up to 
25 /mi (jSoifer et al.lll999l dSN99f ). suggesting them to 
be a major source of the luminosity. Indeed, vigorous 
starburst activity in the nuclear disks was found from 
a concentration of radio supernovae and young super- 
nova remnants within about 50 pc (~0'/2) of each nucleus 
and from radio recombination lines in the nuclear disks 



(ISmith et alj 119981 iRovilos et al. 


20051: lLonsdale et al. 


20061: lAnantharamaiah et al.ll200C 
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20051 ). The presence of an AGN in Arp 220 has not 
been proven nor rejected ba sed on X-ray observations 
because of high obsc uration (jlwasawa et al.|[200ll 120051 : 
IClements et~aLll2002ft . 

Arp 220 appears to represent a short luminous stage 
in merger evolution right before the coalescence of two 
nuclei, because there are luminous merging galaxies hav- 
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ing similar structure. The advanced merger NGC 3256 
(Lg-iooo nm = 10 116 io) has twin gas peaks at two nu- 
clei, possible rotation of the gas around each nucleus, 
and a large (> kpc) g as disk surrounding the central re- 
gion ([Sakamoto et al.il2006l ). NGC 6240 (i 8 -iooo /«n = 
10 119 L Q ) also has two nuclei within 1 kpc on the sky, 
with e ach having a rot a ting stellar cor e and being an 
AGN (jTecza et all 120001 iKomossa et all I2003D . Study- 
ing Arp 220 is particularly important to understand this 
class of objects. 

Accurate spatial distribution of the luminosity should 
tell us much about the hidden energy source(s) in Arp 
220. The dichotomy of the nuclei (with the nuclear disks) 
and the outer disk has been proven useful for this ap- 
proach. S N99| showed using mid-IR imaging and the 
spectral energy distribution (SED) from infrared to mil- 
limeter wavelengths that the two nuclei of ~0'.'3 (~100 
pc) extent have at least 40% of the total luminosity of 
Arp 220. They also pointed out that the outer disk, or a 
kpc-size extended component around the nuclei, could be 
another major component of the luminosity, contributing 
to the SED mainly at around 60-100 fim but too faint in 
continuum to be detected at mid-IR and millimeter wave- 
lengths. Further mid-IR spectroscopy and SED modeling 
suggested, however, that the less-obscured outer disk has 
a moderate starburst with up to 15% of the total lumi- 
nosity, leaving a large m ajority of the lum i nosity to the 
heavily obscured nuclei (ISoifer et ail l2002t ISpoon et al.1 
12004 iGonzalez- Alfonso et alj|2004) The heavily embed- 
de d and compact source of luminosi ty was also suggested 
bv lDudlev fcWynn- Williams! (|1997f ) fr om mid-IR silicate 
absorption. Most of the luminosity in the nuclei is at- 
tributed to starburst by the groups that observed the 
radio supernovae and recombination lines. Still, there 
have been argument s for an AGN(s) (jAalto et al J 120071 ; 
llmanishi et alj|2007l both based on (sub)millimeter line 
ratios favoring X-ray domin ated regions) and even for 
its dominance in luminosity (jHaas et al.ll200ll . based on 
def icit of 7.7 /xm PAH with respect to 850 /jm continuum; 
see lSpoon et a l. 2004 for caveats). 

Recently, IDE07I obtained from their O'.'3-resolution 
observations of 1.3 mm continuum a deconvolved size 
of 0'.'19 x 0'.'13 and a deconvolved peak brightness- 
temperature of 90 K for the western nucleus. They at- 
tributed to the compact source a probable total luminos- 
ity of th e nucle us, 9 x 10 11 Lq in the nucleus-dominated 
model of SN99], and estimated the intrinsic dust temper- 
ature to be 170 K and the 1.3 mm opacity to be ~0.7. 
They argued that the heating source must be an AGN 
on the basis of the high luminosity surface density. Al- 
though this model is possible, the minimum luminosity 
allowed by the 1.3 mm data is 3 x 10 10 L© in the black- 
body limit; it is 1 x 10 11 L Q if the 1.3 mm opacity is 
1.3 as can be estimated from a power-law index o f = 2 
for dust emissivity and the 2.6 mm dust flux in[DE07. 
Thus the majority of the luminosity in the west nucleus 
may come from the starburst in its nuclear disk rather 
than from its continuum core, and the source for the cen- 
tral luminosity may also be massive stars. Our quest for 
the accurate spatial distribution and the source of the 
luminosity in the merger is therefore not yet over. 

In this paper, we report 0'.'2-0'.'3 resolution imaging of 
Arp 220 carried out with the SMA in the submillimeter 



CO (3-2) line and 860 /an continuum. In submillime- 
ter, where dust opacity is higher than in millimeter, we 
can better trace warm dust and better estimate the dust 
temperature and luminosity with less uncertainty from 
opacity correction. This is the first spatially-resolved 
imaging of the merger at a submillimeter wavelength and 
the second submillimeter interferometric observations of 
the galaxy following the original single-baseline exper i- 
ment with the JCMT and CSO ([Wiedner et al.l [20021 ). 
The SMA provided us with comparable or higher reso- 
lutions than previous millimeter observations. Our ob- 
servations are described in [J2] and data reduction in Sj3[ 
The results are presented in []4] for the CO line and in 
SJS] for continuum. Wc analyze the data for the dust 
opacity and temperature, and set a larger lower limit 
to the luminosity in the core of the west nucleus than 
the aforementioned limit from millimeter data by about 
an order of magnitude. The continuum data are also 
compared with the supernova distribution, and the line 
data are analyzed for gas dynamics and ISM properties 
in and around the nuclei. We discuss the properties of 
the hot nuclei and their energy source(s) in f|6j showing 
that some luminosity-related parameters of the west nu- 
cleus are close to the highest for a starburst and that 
submillimeter observations may soon resolve the energy 
source issue. Our finding are summarized in [J7l 

2. SMA OBSERVATIONS 

The Submillimeter Array (SMA) 8 on Mauna Kea, 
Hawaii, was used to observe Arp 220 in 2004-2006. The 
array has eight 6-meter antennas, cry ogenically coole d 
SIS receivers, and a digital corr elator (|Ho et al.1 [20041) . 
The parameters and the log of the observations are in 
Tables [5] and [3[ respectively. 

The redshifted CO(3-2) line was placed in the center 
of the lower side band (LSB) and the continuum data 
were obtained from the LSB line-free channels and from 
the upper side band (USB). The antennas were pointed 
to the center of the merger. Their primary beams are 
approximately Gaussian with a full width half maximum 
(FWHM) of 35" in the LSB. The correlator was config- 
ured for 3.25 MHz resolution and 2 GHz band width in 
each sideband. These correspond to 2.8 km s _1 and 1734 
km s -1 , respectively, for the CO(3-2) line. Data in both 
sidebands were recorded through sideband separation us- 
ing 90° phase switching. 

In each night we observed gain calibrator (s), passband 
calibrator(s), the target galaxy, and, in many cases, the 
flux calibrator Uranus. Either J1613+342 or J1635+381 
or both were observed every 20-30 min during each 
track to calibrate complex gain. Their positions used 
for our observations a gree to 3 mas w ith those from 
VLBI measurements in iMa et ail {l998i) > tying our as- 
trometry to the International Celestial Reference System. 
Passband calibration was made by observing planets and 
bright quasars such as 3C279 and 3C454.3. Pointing was 
checked before and during each track on bright quasars. 

Observations were made on 10 nights but 6 nights were 
lost to poor seeing or various instrumental problems and 
1 night was useable only for flux calibration (track No. 

8 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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5 in Table [3]) . The atmosphere was relatively transpar- 
ent during the four (partially) useable nights; the zenith 
opacity measured at 225 GHz was below 0.1 and hence 
the 345 GHz opacity below 0.4 (see lMatsuo et al.lll998l 
for the conversion). 

The observations in multiple antenna configurations 
gave us projected baselines for the galaxy in the range 
of 10.2 m - 508.1 m after flagging bad data. The range 
corresponds to the fringe spacing of 18" - 0"35 at 345 
GHz. Our longest projected baseline almost matches the 
longest physical baseline of the SMA (508.9 m) because 
Arp 220 transits only 4° from the zenith on Mauna Kea. 

3. DATA REDUCTION 

We reduced our data in three steps. First, the standard 
gain and passband c alibrations were made with MIR (see 
IScoville et al.l ll993. for its original form). Atmospheric 
attenuation was also corrected at this point using the 
system temperature monitored throughout the observa- 
tions. Next, images we re made from the calibrated visi- 
bilities using MIRIAD (|Sault et al.lll995h . Velocity chan- 
nels were binned to 30 km s _1 in this step. We made two 
sets of images, at lower and higher resolution, by using 
different weighting and selection of the visibilities. The 
former has a resolution of 0'.'5, while the latter has a 
resolution of 0'.'33 for the CO line and 0'.'23 for the con- 
tinuum. The high-resolution continuum image used only 
track No. 6 taken in the most extended array configura- 
tion, while the rest of the images used tracks No. 3, 6, 
and 7. Finally, the imag es and data cubes were analyzed 
using the NRAO AIPS (|Bridle fc Greisenlll994D . Details 
of the data reduction are in £)3.1I - £)3.4I 

Throughout this paper, velocities are defined in the 
radio convention, v = c(l — jobs//rcst), with respect to 
the local standard of rest (LSR). All the coordinates 
are in the J2000.0 system. The brightness temperature 
in this paper is the excess temperature above the cos- 
mic background, which is only 3 mJy arcsec -2 at 345 
GHz. Brightness temperatures are calculated with the 
Planck function unless the Rayleigh- Jeans approxima- 
tion is noted or, equivalently, the unit is Krj. 

3.1. Flux Calibration 

The flux scale in our data is based on observations 
of Uranus and its b rightness temperature given by 
iGriffm fc Ortonl ( |1993f ). Flux densities of our gain cali- 
brators were determined by comparison with the planet. 
We separately derived the quasar flux density in the two 
sidebands and averaged the results to obtain the flux den- 
sity at 340-350 GHz, effectively ignoring the small spec- 
tral slope of the q uasars (typically a ~ —0.5 for S„ oc u a ; 
lOear et al.lfl99l . 

The uncertainty of the flux scale of our data was esti- 
mated to be about ±15% by examining relative flux scale 
between our tracks. For this, the flux of Arp 220 in each 
sideband was measured in each track from the common 
baselines in the common ranges of hour angle in order 
to ensure similar u-v sampling. The flux data agreed 
within 15% between the tracks having planet-based flux 
calibration. This method was also used to establish the 
flux scale in track No. 3, which had poor seeing during 
Uranus observations. 

We safely ignored the response pattern of the primary 
beam in all flux and intensity measurements in our data 



because the emission of the merger is compact. Emission 
2" from the center is attenuated by only 1%. 

3.2. Atmospheric Seeing and Astrometry 

We made a test to determine the point spread func- 
tion and to estimate our astrometric accuracy. The data 
from the very-extended configuration (track No. 6) were 
calibrated in the test using only one of the two gain cal- 
ibrators. When calibrating the quasar J1635+381 with 
J1613+342, it had a finite deconvolved size of 0'.'09 in 
FWHM. We attribute this to atmospheric seeing and 
imperfect calibration. (To be precise, the 0'.'09 is not 
exactly the size of the so-called 'seeing disk' because the 
atmospheric decorrelation common to the two quasars is 
calibrated out.) Astrometric accuracy turned out to be 
0'.'03 between the two quasars. Our astrometric error for 
Arp 220 may be larger than this because the galaxy is 
14° and 19° from the quasars while the quasars are only 
6° apart. We used both quasars in the gain calibration 
for the data presented below. 

3.3. Continuum 
3.3.1. Choice of Spectrometer Channels for Continuum 

The continuum data of Arp 220 were obtained by inte- 
grating line-free channels, but the determination of the 
line-free channels needed some care. This is because the 
line emission may have a broad and low-intensity wing, 
and also because there may be weak lines adjacent to the 
target CO line. Recent detections of ~1800 km s~ 1 -wide 
absorption of NH 3 toward Arp 220 (jTakano et al.ll2005h 
and of high- velocit y CO emission in the s imilar luminous 
merger NGC 3256 dSakamoto et al.ll200ll gave us reason 
to be cautious about possible high-velocity gas. 

We adopted as continuum channels those outside of 
the velocity range of 4800-5900 km s _1 after examining 
the CO(2-l) spectrum of Arp 220 shown in Fig. Q] The 
spectrum is from the dual band (690 GH z and 230 GHz) 
SMA observations on 2005 March 2nd (|Matsushita et all 
2008); the array was in the compact configuration. The 
CO(2-l) spectrum has a higher signal-to-noise ratio and 
a wider velocity coverage than our CO (3-2) data, and it 
does not show high-velocity emission outside the above- 
mentioned velocity range. We inferred from this that 
CO (3-2) emission from any high- velocity gas should be 
below our sensitivity limit unless the CO excitation of 
the gas (if it exists) were drastically different from that 
in the main gas component of the merger. 

The continuum data used for the lower resolution (0V5) 
image were from both the lower and upper sidebands and 
were from the same line-free spectrometer channels. This 
choice of the USB channels was a precaution against any 
leakage of the LSB CO line to the USB due to imper- 
fect sideband separation. An error in the SMA software 
had caused a leakage of about —13 dB between sidebands 
until August 2006. For the high-resolution (0'.'23) con- 
tinuum image, we used the USB data of the full 2 GHz 
bandwidth to obtain the highest spatial resolution and 
sensitivity. The CO leakage to the USB should be neg- 
ligible in the data because the line emission was heavily 
resolved out at the resolution. We did not find a differ- 
ence except for sensitivity between the 2 GHz-bandwidth 
map and a test map made only using the band-edge chan- 
nels. We refer to both high and low resolution data as 
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860 /im continuum for simplicity despite their slightly 
different wavelengths. 

Regarding possible contamination of the continuum 
by weak lines other than CO, the large velocity width 
(^1100 km s _1 ) of the galaxy prevented us from notic- 
ing such weak lines overlapping with each other. The 
upper sideband contained a part of the HCO + (4-3) line 
rcdshifted to 350.37 GHz, but its contribution to the con- 
tinuum flux must be small. The expected flux of the 
line is estimated to be 150 Jy km s" 1 from the galaxy's 
HCN(4-3) flux of 260 Jy km s" 1 dWiedner et all 120021 
and the HCO+(1-0)/HCN (1-0) intensity ratio of 0.6 
(jGracia-Carpio et al.ll2006T ) . The continuum flux in the 2 
GHz bandwidth is about 1400 Jy km s" 1 . Thus the con- 
tamination of the USB continuum by the line is expected 
to be below 10% considering that only a part of the line, 
velocities larger than 5200 km s _1 , is in the sideband. 
The actual USB spectrum in Fig. [2] does not convinc- 
ingly show the line, although a hint of excess emission 
may be around 5500 km s _1 . We thus disregard this 
possible line contribution. The upper sideband should 
also contain an edge of the HCN(4-3) line (velocities be- 
low 5000 km s" 1 ), but it can hardly affect the continuum 
flux. 

3.3.2. Continuum Subtraction 

The continuum determined from the LSB line-free 
channels was subtracted from the CO(3-2) line data in 
the same sideband. The subtraction was made in the vis- 
ibilities. The need for this subtraction is not guaranteed 
if there is strong CO emission at the continuum location, 
because the continuum source may be mostly covered by 
molecular gas that is optically thick in CO. An unneces- 
sary subtraction of the continuum would produce an ar- 
tificial depression in the CO maps at the location of each 
continuum peak. In reality, our continuum-subtracted 
CO maps did not show a clear depression. Hence we use 
the continuum-subtracted CO data. The lack of clear CO 
depression at the continuum positions in the continuum- 
subtracted images is indicative of less than complete cov- 
erage of the continuum source by the CO-emitting gas. 
Such partial coverage is consistent with the nuclear disk 
model. The gas coverage at each position would be lim- 
ited to a small velocity range in the model, making the 
velocity-averaged coverage insignificant in the integrated 
intensity maps. 

3.4. Comparison with Single-dish Observations 

The CO(3-2) flux of Arp 220 is (2.2 ± 0.3) x 10 3 Jy 
km s" 1 and (3.2 ± 0.5) x 10 3 Jy km s" 1 in the central 3" 
and 5 ", respectively, in our 0'.'5 resolution data (Fig. [2]). 
The errors are mostly due to flux calibration uncertainty. 
The line flux in the JCMT single-dish observations by 
iWiedner et all (pOOl) is (2.98 ±0.6) x 10 3 Jy km s" 1 for 
the central 14" and (3.7 ± 0.7) x 10 3 Jy km s" 1 for the 
central 22" 9 . Thus the SMA observations detected most 

9 Other singlo-dish CO (3— 2) fluxes in the literature are in 
the range of (3.0 ±1-5) X 10 3 Jy km s" 1 dGerin fc Phillir>sl[T998l : 



of the line emission in the central few kpc of the merger 
though some of the extended emission was likely resolved 
out. The SMA line profile in Fig. [2] matches the JCMT 
spectrum taken at the galaxy center, both being double- 
peaked and having a higher intensity in the blueshiftcd 
peak as in CO(2-l). The CO(3-2) flux detected in our 
0"33 resolution data is (1.6 ± 0.2) x 10 3 Jy km s" 1 As 
expected the high resolution map resolved out more flux. 

The 345 GHz continuum in the central 3" of Arp 220 
is 0.78 ± 0.12 Jy and 0.72 ± 0.11 Jy in our 0"5 and 0723 
resolution images, respectively. Single-dish observations 
of the galaxy at 850 /zm give a flux of about 0.8 Jy with 
a t ypical error of 20% (0.75 Jv ISeiffert et alJl2007t 0.74 
Jv, iKlaas et all 120011: .83 Jy, iDunne et alj |2000[ 0.79 
Jv. lLisenfeld et alj|2000l: 0.83 JvjRigopoulou et alj|1996t 



IMauersberger et a l. 1999; Yao et al.l l2003: Nara yanan et al.1120051) 
The scatter among the observations may be due to, in addition to 
the difficulty in submillimeter calibration, the broad line (~1100 
km s —1 = 1.3 GHz) that is too wide for some spectrometers. 
IWiedner et al. (2002) had a sufficiently large bandwidth with 1.86 
GHz. 



an outlier 0.46 Jv in lAnton et alJl2004f ). Contributions 
of the CO(3-2) line to these bolometer measurements 
are not alwa ys given, but was estimated to be 11% in 
the study bv lKlaas etail (|2001l ). This flux comparison 
indicates that we detected almost all of the continuum 
emission with the SMA. 

4. CO(3-2) LINE IN ARP 220 

Figures [3] and H show our low (075) and high (0733) 
resolution images of the CO (3-2) emission, respectively. 
The two nuclei seen in continuum are marked with plus 
signs. Parameters derived from the data are in Tables H] 
and[5l 

4.1. Spatial Distribution 

The CO(3-2) emission peaks around each of the two 
nuclei (Figs. [3^, and [4^) and has an extended compo- 
nent encompassing the binary nucleus. This distribu- 
tion agrees with those found in CO(1-0) and CO(2-l) 
through previous interferometric observations at l"-073 
resolutions ( *JTJ) - Following SS99J, we call the CO peaks 
associated with the nuclei "nuclear disks" and the ex- 
tended emission an "outer disk" . 

Both nuclear disks appear resolved in our 0733 reso- 
lution CO map with elongation in the north-south di- 
rection. The western nuclear disk appears to have twin 
peaks. The eastern nucleus shows a lopsided CO dis- 
tribution with an extension to the north. This elonga- 
tion agrees with tha t seen in near-IR, where HST images 
(|Scoville et al.lll998f) show two components 'NE' and 'SE' 
in the east side of the merger. These complex distribu- 
tions of CO emission suggest that the nuclear disks have 
non-uniform and non-axisymmetric gas properties. The 
observed distribution of CO(3-2) emission likely reflects 
not only gas distribution but also variation of gas temper- 
ature and line excitation, missing flux, and the blending 
with the outer disk. 

4.2. Velocity Field 

The velocity field of molecular gas in the merger shows 
an overall northeast-to-southwest shift (with NE red- 
shifted) in the outer disk and rapid velocity shifts across 
each nucleus. The apparent velocity gradients across 
the two nuclei are not aligned (Figs. [3Jd and|4})). The 
east nuclear disk has a NE-SW gradient similar to that 
of the outer disk while the west nuclear disk has a 
west-to-east gradient with western side being redshifted. 
The observed velocity structure is consistent with those 
found in previous observations and lead to the model 
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of t wo count er-rotating disks (CO, SS99I; H I absorp- 
tion, IMFP01I: near-IR CO and H 2 , iG^nzelelaLl \Mm 
H53a iRodriguez-Rico et al] I2005T) . The steep velocity 
gradient across the eastern nucleus is also seen in OH 
maser emission (|Rovilos et al.M2003f >. The overall rota- 
tion of the outer disk is consistent with that in previous 
lower-resolution observat ions of molecular gas_jSYB97, 
Downes fc Solomon]|1998h and ionized gas (|Arribas et al.1 
2001H . 

The width of the CO line is clearly peaked at the two 
nuclei (Fig. [3]:) as expected from the steep velocity gra- 
dient across each nucleus. The line profiles at the nuclei 
(Fig. [5J are either flat-topped or doub le-peaked. This 
was first noticed in CO (2-1) bv lSYB9"7l and led them to 
infer the presence of a rotating disk around each nucleus. 
The intensity dip near the line center may be also due to 
absorption by cooler CO near the surface pE07f ). The 
line width at 20% of the maximum intensity in the pro- 
files is 540 km s _1 and 600 km s _1 at the east and west 
nucleus, respectively in our 0'.'5 resolution data. 

The position-velocity (PV) diagrams across each nu- 
cleus confirm the large line widths around, and large ve- 
locity shifts across, the nuclei (Fig. [6|). The line velocity 
changes by several 100 km s _1 within a few 100 pc from 
one side of a nucleus to the other side. In the innermost 
part of the PV diagrams where the apparent velocity 
gradient is nearly linear, we measured a gradient of 2 
km s" 1 pc -1 across the E nucleus and 6 km s" 1 pc -1 
across the W nucleus. The PV diagram of the west nu- 
cleus shows a very large velocity shift in a very small 
region around the center, and he nce the rotation curve 
there may not be close to linear. DE07 suggested that 
the velocity is larger closer to the center in the west nu- 
cleus in their 0'.'3 CO(2-l) data. Although a hint of a 
similar trend is seem in our PV diagram (of 0'.'28 reso- 
lution along the cut) in its redshifted side (offset = 0"- 
0'.'3; 5350-5550 km s _1 ), we conservatively use only the 
velocity shift of 500 km s" 1 in 0'.'22 = 80 pc without 
specifying the shape of the rotation curve. The (appar- 
en t) gra dients that we obtained are consistent with those 
in ISS99I but are larger than the ~0.3-l km s~ 4 pc" 1 in 
previous studies at similar spatial resolutions (jMFPOlt 
iRovilos et all 120031 : iRodri'guez-Rico et al.ll2005h . Those 
lower gradients are probably because they were measured 
in the mean velocity maps, in which gradients are falsely 
reduced through beam smearing. Even our own mean- 
velocity map (Fig. |4] b) gives the apparent gradient of 
about 1 km s" 1 pc -1 for both nuclei, in contradiction 
with the ~500 km s _1 line widths within 0'.'5 (=180 pc; 
Fig. [5j). Also possible, though less likely, is that CO 
and all the other lines in previous studies (H I, OH, and 
H53a) trace different parts of the disks. 

4.3. Brightness Temperature, Column Density 

The peak brightness temperature of the CO line is 
about 50 K. The highest temperatures are seen in the 
vicinity of each nucleus in the high-resolution data (Fig. 
St)- The peak Tb is slightly lower across each nuclear 
disk along the isovelocity contours, seen more clearly in 
the C/33 data though also visible in the 0'.'5 map. These 
temperature troughs are most likely due to a lowered 
filling factor of the molecular gas in each velocity chan- 
nel in the region of largest velocity gradient rather than 
due to an actual decline of the physical gas temperature, 



because the directions of the isovelocity contours are de- 
termined by our viewing angle. Thus the data suggest 
that each nuclear disk has a gas temperature of ~50 Krj 
when averaged over 120 pc (0'.'33) without a correction 
for an area filling factor. The apparent high-temperature 
region is larger around the west nucleus and has a radius 
of about 100 pc at the 40 Krj contour. The gas kinetic 
temperature may be higher if the beam filling factor is 
less than unity or the line is not optically thick. The gas 
in the outer disk shows lower brightness temperatures 
than in the nuclear disks, but it is still mostly higher 
than 10 Krj in the central kpc. The temperature gra- 
dient is probably due to lower filling factor and lower 
physical temperature in the outskirts. 

Fig. [5] shows that the CO line profile and brightness 
temperature are similar in the J=3-2 and J=2-l transi- 
tions at each nucleus. To more accurately measure the 
line ratio, we made matching datasets having the same 
minimum baseline (41 kA) and resolution (0'.'57) from 
our S MA observations and the OVRO observations of 
ISS991 The ratio of integrated brightness temperatures, 
i?( 3 _2)/(2-i), is 1-1 and 1.0 at the east and west nucleus, 
respectively, with an 18% error from flux calibration un- 
certainties. The ratio of unity is expected for optically 
thick emission from warm and thermalized CO. 

The peak CO(3-2) integrated intensity in our data 
is as large as ~1 x 10 4 Krj km s" 1 at the two 
nuclear disks. The column density of H 2 gas is 
Eh 2 = 3 x 10 24 m H2 cm- 2 = 5 x 10 4 M Q pc~ 2 for 
a Calactic conversion factor of Ah 2 //co = 3 x 10 20 
cm~ 2 (Krj km s -1 ) -1 , the applicability of which to 
ULIRGs is not guaranteed. A lower limit to the mass 
surface density of H 2 can be obtained assuming that 
the CO emission is optically thin. For an abundance 
of [CO/H 2 ] = 10" 4 and 50 K LTE excitation, the limit 
is S H2 > 5 x 10 22 m H2 cm" 2 = 7 x 10 2 M Q pc" 2 . The 
excitation temperature of 50 K is a lower limit set by 
the observed brightness temperature; the higher the ex- 
citation temperature, the higher the lower limit of Sh 2 • 
We set a more stringent lower limit later using the CO 
equivalent width ( §6. 2ft . The gas surface density includ- 
ing helium can be obtained by multiplying Ejj 2 by 1.36. 

5. 860 /mi CONTINUUM IN ARP 220 

Figs. [3]and[4]show our 860 /im continuum maps of Arp 
220 at low (0'.'5) and high (0'.'23) resolutions, respectively. 
The continuum has two peaks at the two nuclei and is 
much more compact than the CO line emission. 

5.1. Positions 

The positions of the nuclei in the continuum, listed 
in Table were measured in the two data sets through 
Gaussian fitting and were averaged for each nucleus. The 
positions from the low- and high-resolution maps agree 
to within 0'.'02 for the west nucleus and to 0'.'08 for the 
weaker (and likely more extended) east nucleus. The 
average p ositio ns agree within 0'.'05 with the 1.3 mm po- 
sitions in SS99. The two nuclei are separated by 0'.'98 at 
a position angle of 101°. 

The 860 /im as well as 1.3 mm continuum positions 
agree well with centimeter VLBI observations in terms 
of the relative position between the two nuclei, but show 
a slight (~0'.'l) offset most likely due to astrometric er- 
rors. The median coordinates of the cluster of supernova- 
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related compact (< 1 pc = 3 mas) radio sources around 
each nucleus are a = 15 h 34 m 57!2224, S = +23°30'll'/492 
for the 31 sources around W and a = 15 h 34 m 57!2910, 
S = +23°30Tl'/308 for the 21 sourc es around E These 
are calculated from the so urce lists in lParra et al.1 ()2007l ) 
and lLonsdale et afl (I2006D and are in the coordinate sys- 
tem used~Tn~^a^^e^U |2003) who reported that they 
fixed an error in previous observations. The offset of E 
from W, (Act, AS) = (0?94, -0"18), agrees well with that 
from the SMA. However, both nuclei show an offset of 
O'.'IO ± O'/Ol in almost the same direction between the 
SMA and the centimeter VLBI positions. The offset is 
O'.'ll ± Of! 02 between the |SS99j and the VLBI positions. 
Attributing these offsets to astrometric errors, in which 
ours may be significant (see £)3.2[) . we shift maps to cancel 
the offset when we later compare the spatial distribution 
of the centimeter sources and submillimeter emission. 

5.2. Spatial Distribution 

The two nuclei account for about 75% of the total 860 
/im continuum in the central 3". The rest is in the outer 
disk but its surface brightness is too low for us to con- 
vincingly discern its structure. The flux ratio between 
the two nuclei is about 1 : 2 between E and W. 

Each of the two continuum nuclei is partially resolved 
in our data and has a deconvolved size of ~0'.'15-0'/3 
(Table |4|) . The beam-deconvolved sizes of the nuclei in 
our high-resolution map are larger than the finite de- 
convolved size of a point source, 0'.'09, caused by at- 
mospheric seeing and calibration errors ( £ 13 . 2|) . Thus, 
the submillimeter continuum of each nucleus is intrin- 
sically extended to a few tenths of an arcsecond. The 
deconvolved size of the w est nuc leus from our 0'.'23 map 
agrees with the value that lDE07l obtained from their (X'3- 
resolution map of 1.3 mm continuum. 

The east nucleus has a larger deconvolved size than 
the west and has a major axis roughly in the northeast- 
southwest direction in both of our continuum maps. In- 
terestingly, an elongation similar in size and direction 
is seen and noted in centimeter continuum observations 
(iBaan et all 119871 lNorr"iill988t iBaan fc Haschickl 119951 : 
MFP0¥ lRovilos et all ^OOlTrmd-IR imaging (|SN99h. 
and is also noticeable in the distribution of the radio su- 
pernovae and young sup ernova remnants ([Lonsdale et al.l 
l2006HParra et al J 120071 ) . The west nucleus, on the other 
hand, does not show the east-west elongation seen in the 
distribution of the centimeter compact sources. We dis- 
cuss this in §6.31 

Our imaging results are mostly consistent with the re- 
sults of the previous s ingle-baseline interfero metry at the 
same wavelength by IWiedner et al.1 (|2002f ) . Analyzing 
the visibilities from the 164 m JCMT-CSO baseline, they 
found that the 342 GHz continuum could be modeled as 
a binary system with an east- west separation of ~1" and 
a flux ratio of 1 : 1.7, and that the CO(3-2) line was 
more extended than the continuum. These agree with 
our observations. The total continuum flux attributed 
to the binary nuclei is 0.40 ± 0.09 Jy in the JCMT-CSO 
observations and 0.55 ± 0.08 Jy in our data. Although 
the former makes more room for extended emission when 
compared to the total continuum flux density of about 
0.8 Jy measured with single-dish telescopes, the JCMT- 
CSO visibility data do not hint at such extended contin- 
uum. Namely, the continuum visibilities did not show 



larger amplitudes when the baseline was most foreshort- 
ened and had a frin ge spacing as large as 5" (Fig. 6 of 
IWiedner et al. I [20021 ). This was in contrast to the CO 
visibilities that clearly showed larger amplitudes on fore- 
shortened baselines to suggest extended emission. Thus 
the JCMT-CSO observations and our SMA observations 
agree that the CO emission has a significant extended 
component while the continuum lacks it. 

5.3. Brightness Temperature and Spectral Index 

The two nuclei have high brightness temperatures in 
the continuum, indicating the presence of warm dust. 
The continuum Ty, calculated using the Planck function 
is 27 K and 54 K for the east and west nuclei, respec- 
tively, in the 0H2Z resolution map. The temperatures are 
higher when the sources are deconvolved. A marginally 
resolved source such as these nuclei can be deconvolved 
equally well as a uniform-brightness disk or a Gaussian, 
and the disk diameter is about 1.6 times the Gaussian 
FWHM for the same source (Appendix [A} ■ The west 
nucleus has a FWHM size of ~0'.'15 = 53 pc and a peak 
brightness temperature of 1.6 x 10 2 K in the Gaussian 
model, and a diameter of O'.^S = 85 pc and a temper- 
ature of 91 K in the disk model. The east nucleus has 
a peak temperature of 52 K and 29 K for Gaussian and 
disk models, respectively. The actual dust temperature 
must be even higher than these if the emission is optically 
thin or blurred by atmospheric seeing. 

The spectral index of the continuum emission between 
230 GHz (1.3 mm) and 345 GHz (860 fim) is estimated to 
be a 2 30-34 5 (E) = 2.6 ± 0.6 and a 23 o-34 5 (W) = 2.7 ± 0.5 
for the east and west nucleus, respectively, by comparing 
our S MA observations with the OVRO 1.3 mm data of 
SS99. The index a, defined as S u cx v a , is calculated 
from peak intensities of each nucleus in 0'.'56 x 0'.'51 res- 
olution images at the two frequencies. The SMA map 
used for this was matched to the OVRO data by elimi- 
nating the short baselines (< 40 kA) and tapering long 
baselines. The errors of the spectral indices include those 
from the noise in the maps and the flux calibration uncer- 
tainties. The indices for dust emission only are larger by 
0.1 than the above values. These are calculated by sub- 
tracting the non-dust emission, which we estimate from 
radio data in the literature to be 10 ± 5 (5 ± 5) mJy at 
both 230 and 345 GHz for the west (east) nucleus. 

The spectral indices we obtained at the nuclei are 
larger than that of a blackbody (a < 2) and smaller 
than the index for the entire merger at these frequen- 
cies (a w 3.5; IScoville et al.lll99"lFlDownes fc Solomon! 
1998). The latter index is usually attributed to opti- 
cally thin emission from dust with an emissivity index of 
(3 ~ 1-2 (for e„ oc v&) and one or more temperature 
components (|Eales et all Il989t iRigopoulou et al.l 119961: 
| Lisenfeld et all 120001: IDunne et all 120001: IDunne fc Ealed 
I2OOII : iKlaas et al.1 120011 : iBlain et al.ll2003| ). Our smaller 
spectral indices toward the nuclei are suggestive of mod- 
erate optical depths (i.e., log(ri mm ) ~ 0) toward the nu- 
clei. We attempt a more quantitative assessment below. 

5.4. Optical Depth and Dust Temperature 

5.4.1. West Nucleus 

We constrain the optical depth and dust temperature 
of the west nucleus using the dust sublimation temper- 



860 /im Imaging of Arp 220 



7 



ature, the spectral index around 1 mm, and the com- 
parison of our submillimeter data with mid-IR observa- 
tions. First, a lower limit to the dust optical depth can 
be set from the sublimation temp erature of dust, ~2000 
K (jSalpeteil 119771 : lPhinnevl[l989h . For Arp 220 W, the 
maximum dust temperature and the lower of the decon- 
volved brightness temperatures at the nucleus limit the 
temperature of 860 jitm emitting dust to 90 K < T 8 60 < 
2000 K and dust optical depth at 860 /im to t 860 > 0.05. 

The next constraint on the optical depth, that Tseo is 
of the order of 1 or less, is set from the comparison of 
1300 /im and 860 /im data, as we saw in the previous 
section. The ratio of Rayleigh-Jeans brightness temper- 
atures between two wavelengths Ai and A2 is expressed 
as 

1 - cxp(-r 2 ) 



R 



A2/A1 



l-cxp[-(A 2 /A 1 )/3r 2 ] t 



(1) 



where r 2 is the optical depth at A 2 . Any emission not 
originating from dust has to be subtracted before cal- 
culating the ratio. This formula using the Rayleigh- 
Jeans approximation is accurate to 3% for the minimum 
dust temperature of 90 K and the two wavelengths we 
consider. The ratio R\^/\ 1 is 1 if the emission is op- 
tically thick in both wavelengths, (Ai/A 2 ) /3 if both are 
thin, and in-between for moderate opacities. The 0'.'5- 
resolution OVRO and SMA data of the west nucleus give 
i? 860 /i3oo(dust) = 1.40±0.30. The ratio is also calculated 
from our 0"23 data and the O'.'3-resolution data of DE07 
to be i?860/ 1300 (dust) = 1.57 ±0.32 for the compact west 
nucleus; here we arbitrarily used 10% for the uncertainty 
of the 1300 /im flux including any error caused by differ- 
ent u-v coverage. The limits on the opacities from these 
ratios and Equation (fT]) are summarized in Table [5] For 
example, the former ratio and (3 — 2 give Tg,eo — 2.3 as 
the most likely value and T860 =(1.2-5.3) for the range 
that corresponds to the ±lcr error of i?860/i300- Extend- 
ing this analysis to longer wavelengths is not easy be- 
cause one has less dust emission and more free-free and 
synchrotron emission in the continuum. 

Further constraints can be set on the optical depth and 
dust temperature by combining our submillimeter data 
with the 25 /im imaging photometry as ISN99I demon- 
strated. The model has a bright hot nucleus behind 
a cool shallow absorber whose absorption is significant 
only in mid-IR and emission is insignificant in both w ave- 
length s. Our analysis deviates from that of ISN99I (and 
DE07) in that the nucleus is allowed to have different ap- 
parent sizes between 25 /im and 860 /im, to be consistent 
with the latest observations. The ratio of 860 and 25 
/zm flux densities from the dust can be written as 

£25 _ exp(-T 25 ,f g ) B(25, T 25 ) ^25 ^ 
5*860 1 - exp(-T 86 o) -8(860, T 860 ) ^860 ' 

where S\ and fl\ are the observed flux density of and 
the solid angle subtended by the nucleus, respectively, at 
wavelength A /im, B is the Planck function, T25 is the 
dust temperature at the 25 /im photosphere of the nu- 
cleus, and r 2 s.f g is the for eground opacity that is likely 
to be ~l-2 fSN99l ; see also lHaas et"atll2001l:[Spoon et all 
12004 ) . The nucleus is almost certainly optically thick at 
25 /im because T860 > 0.05 and r 25 ~ t 8 6o (860/25)^ with 
in the range of 1-2. The dust in the compact core 
seen at 860 /im is assumed to be at a temperature of 



about Tseo- Our SMA data and the 25 /im data of SN99 
give S25/SS60 (dust) = 21 and f2 2 5/f2860 = 3.0-7.3, where 
the latter range corresponds to the possible FWHM size 
of the 25 /im emission of 6» 25 =0'.'25-0'/39 (|SN9lf ). Fur- 
ther assuming that the dust at the 25 /im photosphere 
of the nucleus is heated entirely by the more compact 
860 /im-emitting core, we have T25/T860 = 0.76 and 0.61 
for 625 — 0'.'25 and 0'.'39, respectively, since the temper- 
ature would be inversely proportional to the square root 
of radius in this circumstance. These inputs to Eq. ^ 
give an upper limit to the dust temperature in the nu- 
cleus. Namely, T 860 < 1.4 x 10 2 K and 1.5 x 10 2 K for 
625 =0'.'25 and 0'.'39, respectively, if T25,f g is 1. Limits 
for T 2 5 j fg = 2 are in Table [7] The upper limits would be 
lower if there are additional heating sources between the 
860 /im-emitting core and the 25 /im photosphere. For 
completeness, Table[7]also has an upper limit to the dust 
temperature for the case in which the 25 /an and 860 /im 
emission originates from exactly the same surface and 
hence Q.25 = &860 and T25 = Tseo- Because the upper 
limit temperature to the 860 /im-emitting dust does not 
exceed 1.8 x 10 2 K for the plausible combinations of pa- 
rameters, we obtain a lower limit to the 860 /im opacity 
of the west nucleus, r 8 6o > 0.7, from the 90 K lower limit 
to 860 /im brightness temperature. 

The constraints set above are summarized as 90 K 
< ^860 < 180 K and 0.7 < T860 < 5 for the submillimeter- 
cmitting dust in the west nucleus. Given our crude mod- 
eling, in which a single opacity is assigned to an object 
that is unlikely to be a uniform slab, the latter should 
be taken more loosely to mean that T860 of the dust that 
contributes most to the submillimeter emission is of the 
order of 1. We call attention to the fact that the tem- 
perature will be on the lower side of the 90-180 K range 
if there are additional heating sources such as massive 
stars between the bright 0'.'15 core seen at 860 /im and 
the 0'.'25-0'.'39 photosphere at 25 /im. The 860-/im opac- 
ity will then be on the higher side. The moderate optical 
depth in submillimeter toward the nucleus (about 1 when 
converted to r 8 6o ) has been suggested by sev eral authors 
(|SN99t [Gonzalez- Alfonso et al.l[200l IDE07I ). 

5.4.2. East Nucleus 

The east nucleus has less constraining data but its 860 
/im opacity is also likely to be on the order of 1. Here, 
the nucleus has i? 8 60/i30o(d us t) = 1.34 ±0.36 at 0'.'5 res- 
olution. The formal solution of Eq. {T]) is T860 = 0.8 
and 2.8 for (3—1 and 2, respectively. The east nucleus 
has 5 2 5/5g6o(dust) = 16. It is unresolved at 25 /im with 
0'.'62 resolution, but shows elongation or band-dep endent 
positional shift in shorter infrared wavelengths (|SN99l . 
iScoville et al.lfl998h . The foreground extinction to the 
east n ucleus is 50% larger than to the west nucleus 
((SN99T) . This may be the cause of the lower flux ratio 
and color temperature than for the west nucleus. If so, 
the actual dust temperature does not need to be lower 
than in W. For the simplest assumption of O25 = ^ 8 60 
and T25 = i860 j we obtain temperature upper limits of 
(0.9, 1.2, 1.6) xlO 2 K for T 25 ,f g = 0, 1, and 2, respectively. 

5.4.3. Constraint from X-ray Observations 

The optical depth of T860 ~ 1 toward the nuclei corre- 
sponds to a column density of 7Vh 2 ~ 4 x 10 25 cm -2 , and 
a hydrogen mass surface density of ^6 x 10 5 M pc~ 2 
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if the Galactic dust emission law of lHildebrandl (fl983h 
applies. (Our expression of column density in Ah 2 1S f° r 
convenience and not specifying the phase of the hydro- 
gen.) The column density of 10 25 cm' 2 can block the X- 
rays from an AGN that is energetically significant to Arp 
220 and hence coul d explain the non-detection of such 
X-ray s (|Riekel 119881: iDermer et al.lll997t llwasawa et all 
1200 it) . This suggests another possible constraint on rg6o- 
Namely, if T^sjg is larger than the values we used then 
Tsgo could be higher and hence rs6o lower. The bolomet- 
ric luminosity of the brighter west nucleus would then 
exceed 10 12 L@. For example, it would be at least 10 130 
Lq and 10 13 - 7 Lq for r 860 = 0.5 and 0.3, respectively 
f £|5.5[) . An AGN would be needed to explain, among 
other things, a luminosity to mass ratio that is too high 
for a star burst ( §6.3.3j) . A luminous AGN and too low 
T86o would however contradict the X-ray constraint men- 
tioned above, provided that the dust covering the AGN 
emits the bulk of the submillimcter emission. Thus the 
submillimeter opacity is not likely to be much smaller 
than 1, although a precise lower limit is difficult to set 
this way with poorly known dust properties, dust to gas 
ratio, and column density distribution in the nucleus. 

5.5. Luminosity 

The bolometric luminosity of each nucleus can be cal- 
culated from its temperature and size. The luminosity of 
a blackbody sphere of diameter d and temperature T s is 



L, = 7rd 2 aTf 



(3) 



where a is the Stefan-Boltzmann constant. For com- 
parison, the luminosity of a blackbody whose surface 
brightness profile on the sky is a Gaussian of FWHM 
size #fwhm and peak temperature T g is 
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(4) 



if the power we receive is about the average of direction- 
dependent power from the source. As noted earlier, a 
marginally resolved source can be modeled both as a 
uniform-brightness disk and a Gaussian, and the disk 
diameter and the Gaussian FWHM for the same source 
have the relation d « 1-6#fwhm (Appendix |A"|) . The 
(peak) temperatures in these two models have the fol- 
lowing relation 



1 

"m~2 
<0.56T„ 



'FWHM 



T„ 



d < ■" 

(for d = 1.6 6»fwhm) 



The bolometric luminosities have the relation 
L, ps 0.72 L n 



(5) 
(6) 

(7) 



for the same submillimcter flux density, or for the same 
flux density at any frequency where the Rayleigh- Jeans 
approximation applies. In addition, a thin disk of the 
same diameter d has a total luminosity Ld ps 0.5L s if 
we observe it from pole-on, and Ld > L s when edge-on. 
Note that the luminosity estimated from the submillime- 
ter flux density (S v ) strongly depends on the size (R) 
and temperature (T) of the source of any profile as 



L oc R 2 T 4 oc S*iT 8 oc S„T 3 



(8) 



Hence, great care is needed in size and temperature de- 
termination for the luminosity estimate. 

Our 860 /im data give 3.1 x 10 11 Lq and 2.2 x 10 11 Lq 
for the total luminosity of the west nucleus in the Gaus- 
sian and sphere models, respectively, under the black- 
body approximation. The luminosity would be larger 
than that for the sphere if the 860 /zm source was a 
disk and nearly edge-on, which has been suggested for 
the weste rn nuclear disk encompass ing the submillimc- 
ter core (|Scoville et al.1 fl998l [SS99). For comparison, 
the blackbody luminosity is (3-4 ) xlO 10 Lq for the 0'.'3- 
resolution 1.3 mm data of lDE07l . 

These luminosities are lower limits to the true bolomet- 
ric luminosity of the nucleus for two reasons. First, the 
size used to estimate the luminosity may be larger than 
the true size because a point source would have a finite 
size of 0'.'09 in our data due to atmospheric seeing and 
various calibration errors ( ^3. 2[) . Removal of this would 
make the source FWHM, peak temperature, and the to- 
tal luminosity 0" 11 = 40 pc, 2.7 xlO 2 K, and 1.4xlO 12 L , 
respectively, for the Gaussian model. We warn however 
that these numbers, in particular the luminosity, have 
large uncertainties and that they should be taken only 
as possible values until supported by higher-resolution 
observations. The second but more important reason 
why the (2-3) xlO 11 L© is a lower limit is that our 860 
/jm brightness temperatures are lower limits to the true 
temperature of the dust because of the moderate opacity 
of the dust emission. The dust temperature would be 1.6 
times higher than the observed brightness temperature 
if 7860 is 1, and that would make the bolometric lumi- 
nosity ^6 times larger to L ~ 10 12 L©. The use of the 
Stefan-Boltzmann law in this calculation is valid because 
most of the luminosity is emitted in mid- to far-IR where 
the dust is opaque 10 . Although the opacity correction in- 
creases the luminosity by only 3% if T86o = 5, it is almost 
certain, as can be inferred from our CO maps, that the 
dust around the nucleus has a lower-column density pe- 
riphery that is transparent at 860 /im but is opaque at 
shorter wavelengths. This luminosity from outside of the 
bright 860 /im core but within the nuclear disk is proba- 
bly going to increase the luminosity of the nucleus to the 
level expected in the current model of luminosity distri- 
bution in the merger system (JO). The extra luminosity 
of the periphery would be either reemission of the radia- 
tion from the hot core or the luminosity of the starburst 
in the nuclear disk. 

The east nucleus has a (minimum) luminosity of ^5 x 
10 9 L Q before any seeing and opacity corrections. The 
luminosity of the outer disk is poorly constrained by our 
continuum observations. Our line data, however, suggest 
that it may be as large as a few 10 n .L© as in the cur- 
rently prevailing model described in ^JTJ The outer gas 
disk of ~1 k pc diam eter and several 10 3 M© pc -2 sur- 
face density (|SYB97ft is opaque shortward of ~50 /im and 
will have a total luminosity of a few 10 11 Lq if it is at 40 
K. The high peak brightness temperatures of the CO(3- 
2) line, > 30 K, seen in the central 2" (0.7 kpc) of the 
merger hint at such a warm ISM in the outer disk. In ad- 

10 The luminosity of a gray body with a frequency-dependent 
opacity of t(u) = (yjv\)P (/3 =1-2) is more than 90% of the black- 
body luminosity for the same temperature if hv\/kT < 1, where 
u\ is the frequency at which r is 1. 
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dition to the luminous nuclei, the outer di sk has its own 
young massive clusters for heating (e.g., IScoville et al.1 
1998; lWilson et alJl2006h . 



6. PROPERTIES OF THE NUCLEI 
6.1. Mass and Dynamical Properties 

The rotating disk around each nucleus is the currently 
accepted scheme, although there were alternative models 
attributing the velocity struc ture to velocity dispersion 
or non-circular motion (e.g., iDownes fc Solomonl 119981 : 
lEckart &; Down cs 2001). For the west nucleus, for ex- 
ample, HST/NICMOS images imply an inclined disk of 
absorptio n (i.e., gas and dust) elongated in the east- west 
direction (|Scoville et all 119981) . The major axis and the 
extent of the disk agree with those of the observed gas 
velocity gradient. Radio supernovae and their remnants 
in the west nucleus are in a 160 p c by 70 pc region 
elongated in the eas t -west direction (ISmith et alj 119981 : 
lLonsdale et al.l 120061 : iParra et"aT] l2007f ). also consistent 
with the inclined and starbursting nuclear disk elongated 
in the same direction. IDE07I also adopt a rotating disk. 

The dynamical mass of each nucleus is ~10 9 Mq within 
~100 pc radius (see §4.21 and Tabled]) as was prev iously 
derived from CO( 2-l) and near-IR observations (SS99, 
iGenzel et alj |2001). The estimate should be accurate to 
within a factor of a few, but it is difficult to tell in which 
direction it is likely to err. On one hand, the calculated 
mass underestimates the true mass because no inclina- 
tion correction is applied, though the correction is likely 
small for the highly inclined west nuclear disk. On the 
other hand, the radius at which the rotational velocity 
peaks in each PV diagram is quite possibly overestimated 
because each nuclear disk is only marginally resolved in 
space; the dynamical mass will be overestimated by this 
error. We can not tell which of the two sources of bias is 
larger, and hence the calculated mass is neither a lower 
limit nor an upper limit. Still, the current data do not 
support a small mass of ^2 x 10 8 / sin 2 i Mq within a bout 
100 pc from the west nucleus suggested bv lMFPOU . The 
small velocity gradient seen in H I absorption that lead to 
this mass estimate is incompatible with the large CO ve- 
locity range seen in our PV diagram (i j4.2[) . Their model 
of merger configuration and history are partly based on 
the observation that the west nucleus is an order of mag- 
nitude lighter than the east nucleus, but that disparity 
is not supported by our data. If there is any disparity, 
the west nucleus appears to be denser than the east nu- 
cleus as suggested from the larger velocity gradient of 
the former (see below), but that could be due to a lower 
inclination of the eastern disk. 

The mean mass density calculated from the V/R and 
the following formula, which does not assume a linearly- 
rising or any specific rotation curve, is 2 x 10 2 Mq pc~ 3 
for r < 110 pc of the east nucleus and 2 x 10 3 Mq pc~ 3 
for r < 40 pc of the west nucleus. 

„-g-V)=55(- V >* )' (9) 
^Mopc" 3 / V kms /p<v 

Elliptical galaxies a nd bulges have sim ilar mass densities 
at the same scales (|Lauer et al.l 119951 ) an d our Galaxy 
has o (r = 12 pc) = 4 x 10 3 M pc -3 (jGenzel et all 
1996). The mass densities that we derived are lower lim- 
its, because both the underestimate of the rotational ve- 
locity V due to the lack of inclination correction and the 



overestimate of the radius R due to insufficient spatial 
resolution make the ratio V/ R and hence the mean mass 
density smaller. For the same reason, we obtain upper 
limits of 1-3 Myr to the dynamical time-scales (= rota- 
tion periods) of the nuclear disks from the velocity-radius 
pairs. 

6.2. ISM Properties Suggested from Low CO(3-2) 
Equivalent Width 

In addition to the high brightness temperature, our ob- 
servations show that the nuclear disks have a much lower 
line-to-continuum ratio than the oute r disk in the submil- 
limeter, as was the case at 1.3 mm (SS99). Specifically, 
about 75% of the total (i.e., single-dish) 860 //m contin- 
uum comes from the nuclear disks while only about 25% 
of the total CO(3-2) emission is from the nuclear disks. 
The rest of the emission is from the outer disk. To put it 
another way, the equivalent width for the CO(3-2) line 
is ~1 x 10 3 km s _1 in the nuclear disks, while that of the 
outer disk is several times larger, ^6 x 10 3 km s _1 . 

The larger equivalent width of the outer disk is in line 
with the observations of other IR bright galaxies and 
the small equivalent width (i.e., low line-to-continuum 
ratio) in the nuclear disks is exceptional. For example, 
ISeaauist et ail (|2004f ) showed that CO(3-2) emission typ- 
ically comprises 25% of 850 /jm flux detected within a 
~3 x 10 4 km s~ 4 bandwidth at the central kiloparsecs 
of relatively IR-luminous galaxies (Io^,(Lib./ Lq) ps 10- 
12) in the SLUGS survey (jDunne et aLll2000t lYao et all 
12003( 1. The fraction of 25% means a CO (3-2) equiva- 
lent width of 1 x 10 4 km s~ 4 . Similarly, the equivalent 
widths in the central 15" (250 pc) of the starburst galax- 
ies NGC 253 and M82 are (1.5 ± 0.3) x 10 4 km s" 1 and 
(l.l±0.2)xl0 4 km s -1 , respectively, a ccording to the 345 
GHz observations of these gal axies (jAlton et all 119991 : 
Israel et al.lll995l: lTilanuslll991f ) and the prescription by 



Seaquist et all (|2004h to remove the CO(3-2) contribu- 



tion from the 850 /im bolometer flux. 

A simple model tells us that the low equivalent width is 
basically due to either high temperature or high column 
density per velocity. Figure [7] shows LTE calculations 
for CO(3-2) optical depth, excess brightness tempera- 
ture, and equivalent width. The model assumes that gas 
and dust in the model region have the same uniform tem- 
perature, gas-to-dust mass ratio of 100, CO abundance 
of [CO/H2]=10 -4 , and dust mass opacity coefficient at 
the CO(3-2) frequency of k = 1 (g cm" 2 )" 1 , and also as- 
sumes that the dust emission is optically thin. Appendix 
[B]gives model formulae and the reason why the moderate 
dust optical depth of the nuclear disks does not affect the 
derived gas column density much. Figure 7 shows that 
there are two possible sets of gas conditions for each set 
of equivalent width and brightness temperature. One has 
optically thin CO emission from warmer gas with lower 
column density per velocity, and the other has optically 
thick CO emission from cooler gas with higher column 
density per velocity. 

The optically thin solution is unlikely for the nuclear 
region of Arp 220 because it needs a high temperature of 
about 800 K for the nuclear disks. The required tem- 
perature would be even higher for plausible non-LTE 
cases in which levels much higher than J=3 are less popu- 
lated than in LTE because of their high critical densities. 
Such high temperatures are incompatible with the upper 
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limit set on the dust temperature in i )5.4l Moreover, 
the CO(3-2) to (2-1) intensity ratio i?( 3 _ 2 )/(2-i) is close 
to 1 at the two nuclei, consistent with optically thick 
emission. The ratio would be | if both lines were op- 
tically thin. Accepting that the CO (3-2) emission from 
the nuclear disks is optically thick, the bulk of the CO 
emission from the outer disk must also be optically thick 
because the intensity ratios b etween low-J CO lines are 
about 1 for the entir e galaxy fMaucrs berger "eTaT][l99l 
IWiedner et al.|[2002T ). Although the optically-thin solu- 
tion of 800 K gas is unlikely, it provides a lower limit 
to the mean gas surface density in the nuclear disks, 
Sh 2 > 4 x 10 23 m H2 cm~ 2 = 6 x 10 3 M© pc~ 2 , because 
this solution has less gas surface density than the opti- 
cally thick solution. This is a more stringent lower limit 
to the column density than the one in thanks to the 
additional constraint from equivalent width. 

In the preferred case of optically thick CO emission, 
the low equivalent width is due to a high column density 
per velocity width, iV C o/AV w 10 18 5 cm -2 (km s -1 ) -1 
for the equivalent width of 10 3 km s^ 1 . The H2 col- 
umn density in the nuclear disks would be Eh 2 = 1 x 
10 25 m H2 cm" 2 = 2 x 10 5 M Q pc -2 if the line width in 
the disk is 300 km s -1 , which is half of the total line 
width due to disk rotation and is an upper limit to the 
local line width in the disk. The high column density is 
in order-of-magnitude agreement with what we inferred 
from the dust emission alone in H5A[ as it should be. The 
high column density A/h 2 ~ 10 25 cm~ 2 suggests that the 
volume density in the nuclear disks is high, nn 2 > 10 4 ' 5 
cm -3 for a volume filling factor of 1 and a line-of-sight 
depth of < 100 pc for the disks. The high density 
of the molecular gas is in agreement with the observa- 
tions of molecular lines with high-critical densities (e.g., 
IWiedner et alj|2002t llmanishi et al.ll2007l ) and is also in 
line with models of H II regions and supernova r emnant s 
(jAnantharamaiah et al]|2000t iParra et al.ll2007| ). 

A more elaborate model that also employs optically 
thick CO emission and explains the lower CO equivalent 
width in the nuclear disks is that the ISM has a tem- 
perature gradient in such a way that the gas and dust 
are hotter inside and cooler near the surface of the disk. 
This reduces the CO emission through self-absorption 
and hence reduces the equivalent width of the line, be- 
cause the CO line has higher optical depth than the con- 
tinuum, Tsgo ~ 1. The absorption feature in the CO(2-l) 
line toward the west nucleus, observed bv lDE07l . suggests 
that this is a part of the reason for the low equivalent 
width. It is quite possible that the disks have not only 
the temperature gradient but also a density gradient with 
the center having higher column and volume densities, as 
hinted by the compact dust emission. 

6.3. Energy Source 
6.3.1. Sources in the Nuclear Region 

Figure [8] shows the centimeter radio sources and the 
submillimeter continuum in the nuclear region. Relative 
positions of various centimeter sources in the figure are 
based on observations but the registration of centimeter 
and submillimeter data is made, as noted in §5.11 with 
an assumption that the two nuclear sources seen in both 
wavelengths are the same. Specifically, the registration 
made the submillimeter peak coincide with the median 



centroid of the supernova-related sources (hereafter SNe) 
in the brighter west nucleus. One can alternatively use 
the peak of diffuse centimeter emission for the registra- 
tion, because the offset between the two nuclei in that 
emission, (Aa, AS) = (0"96, -0"15) (jBaan fc Haschickl 
1 1995T ) , also agrees with the offset measured in submillime- 
ter. However, the western peak of the centimeter diffuse 
continuum (green circle in Fig. [8]) is only 40 mas from 
the center of the SNe distribution (red circle in Fig. [8]) 
according to the relative astromet ry of the diffu s e emis - 
sion with respect to the SNe by iRovilos et al.1 |2003). 
Thus our centimeter-submillimeter registration in Fig. 
M is probably accurate to MX' 04. 

Neit her of the two OH m egamasers in the west nu- 
cleus ( Lonsda le et al.l 1 19981 . shown in blue in Fig. [8]) 
coincides with the submillimeter continuum peak, ac- 
cording to the registration made above. Of particu- 
lar interest is t he sou thern maser source called W2 in 
lLonsdale et al.l (|1998l ). It has a large velocity gradient, 
dV/dr w 18.7 km s _1 pc -1 in a ^6 pc diameter region, 
and may trace a hea vily obscured AGN of ~1.7 x 10 7 M© 
(jRovilos et all 120051) . Whether it is an AGN or not, our 
submillimeter map and its registration in Fig. [8] suggest 
that it is unlikely to be the main heating source of the 
submillimcter-emitting dust. We checked if the submil- 
limeter data could be registered with the centimeter data 
using the OH megamasers and found it unlikely Among 
the four OH megamasers in the nuclear region, the pair 
between the southern source in each nucleus has the sepa- 
ration that is closest to the one between the submillime- 
ter peaks. This maser separation, (Aa, AS) — (0'.'89, 
— 0'.'15), is in worse agreement with the submillimeter 
separation (0'.'96, — 0'.'18) than the supernova-based sep- 
aration (0'.'94, — 0'.'18) and the abovementioned separa- 
tion in diffuse centimeter continuum. Moreover, the shift 
needed for the registration is larger if the southwest 
megamaser should coincide with the west submillime- 
ter peak. Furthermore, the elongation in the northeast- 
southwest direction seen both in the submillimeter emis- 
sion and the supernova distribution of the east nucleus 
suggests that these two, rather than the submillimeter 
and OH emission, should be correlated. Thus, although 
the final conclusion should be made from more accu- 
rate astrometric observations, we favor the registration 
in Fig. [5] and the SNe-submillimctcr connection over the 
megamaser-submillimeter connection. The maser emis- 
sion off the west nuclear disk may be collisionally excited 
by a bipolar wind or outflow from the disk. 

6.3.2. Supernovae and Dust Emission 

The spatial distribution of the submillimeter contin- 
uum is expected to follow the surface density of super- 
novae under certain circumstances that include the dom- 
inant energy source being a dust-enshrouded starburst 
(see Appendix lCj). We show the smoothed distribution 
of the SNe in Figure [51 In Fig. [5^ the smoothing ker- 
nel is a 0'.'5 Gaussian that is about the same as the PSF 
of our low-resoluti on da ta (Fig. [3]) and the PSF of the 
mid-IR images of [SN99J. As noted earlier, the eastern 
nucleus has an oval shape elongated in the same NE-SW 
direction in the supernova distribution, 860 fim contin- 
uum (see the deconvolved shape in Tabled]), and mid-IR 
continuum, among others. This agreement is suggestive 
of the dust being heated mainly by the starburst, al- 
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though there may still be room for an accreting black 
hole to hide within or behind the starbursting disk. In 
contrast to the east nucleus, the brighter west nucleus 
shows different shapes in the submillimeter continuum 
and in the supernova distribution (see Figs. [9Jd and[8|). 
The dust continuum is more compact and the starburst 
is much more elongated in the east-west direction at the 
0'.'23 resolution. This bright 860 /mi core has a minimum 
luminosity of 2 x 10 11 Lq, calculated in M 5 . 51 

The reason for the different spatial distribution be- 
tween the submillimeter continuum and SNe must be 
either that many massive stars and their remnants are 
missed in the current VLBI observations at the densest 
region or that there is a dust-heating source other than 
massive stars in the center of the west nuclear disk. The 
former could be because of free-free absorption, radio 
quiet supernovae, younger stellar population at the cen- 
ter of the W nucleus than the one around it and in the 
E nucleus, or more/less dust and gas in the nucleus/disk 
of Arp 220 W than its star formation implies. The non- 
stellar heating source, if it exists, must be an AGN. 

Quantitatively, the west nucleus has three of the four 
new radio sou r ces, i.e., radio supernovae, found by 
lLonsdale et~aTl (|2006f) in the merger over a period of 
one year. T wo of them are in the center of the west 
nucleus, and iParra et alJ (|2007f ) pointed out that more 
core-collapse supernovae are probably left undetected be- 
cause the detected ones are most likely type Iln that 
are radio-luminous and rare. It may therefore be pos- 
sible that there is a stronger concentration of super- 
novae resulting from vigorous star formation at the cen- 
ter of the west nuclear disk than is visible in the existing 
VLBI observations. The starburst that caused the ob- 
served supernovae can provide the entire luminosity of 
the merger even before the correction for non-type Iln 
SNejSmith et alll998t lLonsdale et~alll2006t IParra et all 
l2007f l. although such an estimate involves assumptions 
on the stellar mass function and starburst history, which 
will be difficult to verify. 

6.3.3. Constraints from Luminosity-to-Mass Ratio and 
Luminosity Surface Density 

Among quantitative parameters that could constrain 
the energy source from our data are luminosity-to-mass 
ratio and luminosity surface density. A lower limit to the 
luminosity-to-mass ratio within r = 40 pc of the west 
nucleus is about 4 x 10 2 Lq/Mq according to the pa- 
rameters estimated in §5.51 and £ 16.11 The ratio would be 
as high as 10 3 Lq/Mq if the nucleus has a luminosity of 
10 12 Lq as shown possible depending on its size and opti- 
cal depth. The L/M ratio, which depends on stellar mass 
function, may be as high as ~1 x 10 3 Lq /Mq for a very 
youn g starburst (age < 10 Myr, see, e.g.. lLeitherer et all 
1999). Its more robust upper limit is ~1 x 10 4 Lq/Mq 
for a pure population of 100 M© stars. The high ratios 
in the range of (10 3 -10 4 ) Lq/Mq are observed in super 
star cl usters (SSCs) such as the one in NGC 5253 stud- 
ied bv lTur ncr et al. (200 0]) an d the Quintuplet and other 
clusters in Fig er et all ( 1999h . For comparison, an AGN 
fueled at the Eddington rate has a luminosity to black- 
hole mass ratio of 4 x 10 4 Lq/Mq. If the L/M ratio is 
close to our lower limit, then it is quite possible, as far 
as this parameter is concerned, that the main luminosity 
source of the nucleus is an intense starburst producing 



many super clusters. The L/M ratio for the 10 12 Lq 
case, on the other hand, is close to the largest a star- 
burst can have, considering that the observed dynamical 
mass includes everything in the region, such as young 
and old stars, massive black holes (if any), and the gas 
and dust of high column densities. We can therefore in- 
fer that the west nucleus, if it has a luminosity of ~10 12 
Lq , must have either an energetically significant AGN or 
an extreme starburst that is compact (r < 40 pc) , young 
(age < 10 Myr), and maybe biased toward high mass 
stars. 

The luminosity surface density of the west nucleus is 
10 7 - 6 Lq pc" 2 or lar ger in the central 80 pc. The lower 
limit is already am ong the highest fo r infrared-luminous 
starburst galaxies (jSoifer et all 120011 their Table 4). If 
the west nucleus has a luminosity of 10 12 Lq within 50 
pc diameter, then the parameter reaches 10 87 Lq pc~ 2 . 
Obviously, a luminous AGN would have no problem to 
achieve this level of luminosity surface density. At the 
same time, there are known stellar systems that could 
also achieve the high values. Some of the youngest su- 
per star clusters have core radii of 0.1-1 pc, luminosities 
of ~10 7 -10 9 Lq, and ag es less than 10 Myr (jFiger et all 
[1991 lTurner"etaII[2000h . Taking a luminous SSC of 10 9 
Lq as a yard stick, there would be 200 and 1000 of them, 
respectively, for the (L, d)=(2 x 10 11 Lq, 80 pc) and 
(10 12 Lq, 50 pc) cases mentioned above if the SSCs are 
the main luminosity source. Their mean spacing would 
be 11 pc and 4 pc, respectively. The clusters would prob- 
ably overlap in their outskirts, since half of the ~7000 O 
stars in the most luminous SSC in NG C 5235, used as 
the y ard stick, are in the central 5 pc (|Turner fc Beckl 
2004) . Alternatively, many of the clusters may have al- 
ready been disintegrated to form a single concentration 
of massive stars in the galaxy nucleus, consid ering the 
short dissolution time of SSCs (|Fall et al.l l2005). It could 
also be that a smaller number of heavier cluster(s) was 
formed rather than hundreds of SSCs. This possibility 
cannot be excluded on the grounds of the number density 
of massive stars required to produce the luminosity sur- 
face density of Arp 220 W, since higher number densities 
are observed in SSCs. 

6.3.4. Star Cluster v.s. AGN 

As seen above, the parameters of the luminous west 
nucleus are in or close to the limit explicable by a star 
cluster or clusters. The starburst model may explain 
some of the observations that appear to favor an AGN 
but it certainly has some difficulties. The extreme lu- 
minosity surface density even for a ULIRG is expected 
to result in photodestruction of PAH particl es as was ob- 
serve d around the model SSC in NGC 5253 (|Beirao et al.l 
2006). This may be part of the reason for the merger's 
low PAH emission to 850 /im c ontinuum ratio th at was 
taken as evidence for an AGN (Haas et~aT][200l . The 
excess dust emission at the center of the west nucleus 
compared to the number of supernovae and young su- 
pernova remnants is reminiscent of the abovementioned 
SSC in NGC 5253 whose predominantly thermal radio 
spectrum is attr ibuted to a very young starburst with 
fe w supernova e (|Beck et all 119960 . On the other hand, 
as lHaas et "all (|2001[ ) emphasized, a cluster of stars (or 
SSCs) with an extent of a few 10 pc or larger would 
be more difficult to enshroud with dust than an accret- 
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ing black hole. The cluster model may therefore be in- 
compatible with the deeply ens hrouded but energet ically 
dominant component used by ISpoon et alj (|2004fl and 
iDudlev fc Wvnn-Williamsl (fl997T to model the mid-IR 
spectrum. Another difficulty in the cluster model is that 
it needs to form hundreds of SSCs or equivalent high- 
mass stars at high volume density in the short time of 
< 10 Myr. The short time scale is necessary because the 
luminosity-to-mass ratio of a cluster declines once mas- 
sive stars start to die and also because gas and dust in the 
nucleus would be dispersed by the energy injection from 
the starburst. Alt hough SSCs tend to form in groups 
(jZhang et al.ll200"ih . we do not yet know the mechanism 
for such a compact and intense burst of SSC formation 
or even whether it is possible. 

A massive black hole of ~10 8 M Q can generate (0.2- 
l)xl0 12 Lq via mass accretion at a sub-Eddington rate, 
as it is believed to do in quasars. The above-mentioned 
constraints on mass-to-luminosity ratio and luminosity 
surface density of the west nucleus can be easily met with 
a buried quasar. The weaker concentration of supernovae 
and supernova remnants toward the center of W than the 
concentration of the submillimeter continuum and bolo- 
metric luminosity can be explained if a large part of the 
luminosity is due to an AGN. There is enough material, 
Ah 2 ~ 10 25 cm -2 , to block the X-rays, and it is easier to 
enshroud an AGN than a million O stars. An AGN with 
a high-covering factor could explai n the large equivalen t 
width of the Fe Ka line of Arp 220 (jlwasawa et al.ll2005[ K 
Regarding the formation mechanism, it is not surprising 
that the remnant nucleus of one or both of the merged 
galaxies has a massive black hole, since most if not all 
galaxies are thought to harbor a massive black hole at 
the center. The massive black hole in Arp 220 W would 
be at the center of the nearly edge-on nuclear disk of gas 
and dust, having abundant gas in its vicinity to accrete, 
and might have been significantly growing through the 
increased supply of fuel. The number of supernovae and 
the amount of free-free emission that appear to suggest a 
starburst large enough to provide the entire luminosity of 
Arp 220 is a problem for the luminous- AGN model. Our 
observations suggest that neither of the two OH mega- 
masers about 50 pc off the luminosity center marks the 
main luminosity source even if they are associated with 
AGNs. 

Taking everything above into account, a luminous 
AGN at t he cent er of the west nucleus, proposed most re- 
cently by DE0 7|, seems to have fewer difficulties than the 
SSC-cluster model if the luminosity of the bright 860 fim 
core is significantly larger than the lower limit we set. If, 
on the other hand, the luminosity is close to the observa- 
tional lower limit, then the majority of the luminosity of 
the west nucleus is emitted from the nuclear disk rather 
than the core. In this case the heating source for the 
disk, as well as for the core, can be a starburst. In other 
words, the spatial distribution of luminosity within the 
west nucleus is fo und m ore centrally peaked than known 
before (including lDE07h , but the parameters of the con- 
centration obtained so far are not enough to conclude 
the existence of an energetically-dominant AGN. Really 
robust evidence is needed here considering the number 
of previous studies that have concluded one way or the 
other about this issue. 



6.4. Prospects 

Further observations will let us set tighter constraints 
on the energy sources and their spatial distributions. For 
example, higher resolution observations at around 860 
fim will better determine the luminosity-to-mass ratio. 
Such observations may soon become possible through the 
ongoing eSMA project that combines the SMA with the 
neighboring JCMT and CSO telescope. Observations 
at even shorter submillimeter wavelengths, where dust 
opacity is higher, will better constrain the dust tempera- 
tures in the nuclear and outer disks and will thereby de- 
termine their relative contributions to the total luminos- 
ity. They also let us search for a compact and bright core 
in the east nucleus. The luminosity estimate from dust 
thermal emission in the submillimeter wavelengths com- 
plements other methods, because it bypasses the number 
of ionizing photons or supernovae to estimate luminosity 
and also because it is less hindered by foreground extinc- 
tion than infrared-based methods. In addition, various 
molecular lines in millimeter and submillimeter will tell 
us the physical properties of the ISM and alert us to 
any anomaly in the nuclei. Ongoing projects in these di- 
rections, further modeling of the structure and radiation 
of the ISM, and the eventual arrival of ALMA will tell 
us the heating sources not only in Arp 220 but also in 
ULIRGs in general. As an example, high-resolution pho- 
tometry at a ~10 pc resolution across the ALMA observ- 
ing bands will allow us to locate a Compton-thick AGN 
of 10 12 L & as a source of a few 100 K, even if it is hidden 
in other wavelengths behind layers of absorbers. The dy- 
namical mass estimated from molecular lines may reveal 
that its luminosity to mass ratio is too large for a star- 
burst. Finding many such sources may prove the long- 
hypothesized formation of quasars in lum inous mergers 
([Sanders et alJll988t I Hopkins et al .1120061 and references 
therein) . 

7. SUMMARY 

The ultraluminous merging galaxy Arp 220 has been 
imaged using the Submillimeter Array (SMA) at Of! 2- 
0'.'3 (-100 pc) resolutions in the CO(3-2) line and 860 
jjxa continuum. We analyzed structure and properties 
of the submillimeter emission with particular attention 
to the spatial distribution of bolometric luminosity. The 
luminosity distribution allows us to constrain the energy 
source of Arp 220. Our main results are as follows. 

1. Both the CO(3-2) line and the submillimeter con- 
tinuum peak at or around the two nuclei of the merger. 
The continuum emission is mostly from compact (a few 
O'.'l) regions at the two nuclei. The CO gas is more ex- 
tended than the continuum at the nuclei, and it also has 
a component with a total extent of ~3"(1 kpc) that en- 
compasses the double nucleus. The overall distribution 
of CO and dust emission is in agreement with the obser- 
vations at lower frequencies. 

2. The peak brightness temperature exceeds 50 K in 
both the CO (3-2) line and the submillimeter continuum 
emission, revealing the presence of warm gas and dust 
in the center of the merger. The west nucleus has the 
strongest continuum emission that is consistent with a 
Gaussian of FWHM O'.'l 5 (50 pc) and peak brightness 
temperature 1.6 x 10 2 K or a uniform-brightness disk of 
diameter 0'.'23 (80 pc) and a brightness temperature of 
0.9 x 10 2 K. 
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3. The dust opacity at both nuclei is on the order of 
unity at 860 /im, according to the analysis of flux ratios to 
1.3 mm and to 25 /im. Both the 860 /mi opacity and low 
equivalent widths of the CO (3-2) line at the two nuclei 
(~10 3 km s^" 1 ) suggest a high column density of gas and 
dust toward the nuclei, of the order of JVh 2 ~ 10 25 cm~ 2 . 
The equivalent-width analysis and the CO (3-2) to CO (2- 
1) intensity ratio close to 1 suggest that the CO(3-2) line 
is optically thick and the molecule is thermalized at least 
up to the J=3 level in the nuclear disks. The gas there 
must be also dense, nn 2 > 10 4 5 cm~ 3 . 

4. The bolometric luminosity of the bright 50-80 pc 
core within the west nucleus is at least 2 x 10 11 Lq and 
may be 10 12 Lq or higher depending on the opacity and 
structure of the source. Our blackbody lower limit is 
an order of magnitude larger than that from existing 1.3 
mm observations. The east nucleus has a lower limit to 
the bolometric luminosity of 5 x 10 9 Lq. It also may be 
more luminous if its 860 /im emission is optically thin. 

5. The sharp velocity shift across each nucleus is con- 
firmed in the CO (3-2) velocity field, and so is the overall 
rotation of the larger (~kpc) gas disk. Each nucleus has 
a line-width of ^500 km s _1 and a dynamical mass of 
~1 x 10 9 Mq within a radius of ^100 pc if the gas mo- 
tion is due to rotation, as seems most likely. The mean 
mass densities in the nuclei are on the order of 10 2 -10 3 
M & pc -3 , and are comparable to those of galaxy centers 
at the same scale. The overall merger system has a full 
CO line width of ~1100 km s _1 and does not show ex- 
tremely high velocity CO at our sensitivity. The overall 
structure and kinematics are consistent with what has 
been proposed, a binary system of counter-rotating nu- 
clear disks encompassed in a larger disk. 

6. Morphological agreement between dust continuum 
and surface density of supernova-related VLBI sources 
in the east nucleus is suggestive that the dust is mainly 
heated by a starburst. In contrast, the west nucleus has 
more compact dust emission than the supernova distribu- 
tion. The compact hot core has a luminosity-to-mass ra- 
tio of > 4 x 10 2 Lq/Mq and a luminosity surface density 
of > 10 7 6 Lq pc~ 2 in the central 80 pc, with a possibil- 
ity of actually having ~5 times larger values than these 
lower limits. These parameters, at least the lower limits, 
are similar to those in young super star clusters and are 
close to the highest for a starburst. Hence we can not yet 
rule out an extreme starburst equivalent to hundreds of 
young SSCs. A buried luminous AGN accompanied by a 
starburst can also account for the luminosity and other 
observed parameters of the west nucleus, and it becomes 
more plausible in higher luminosity. Further observa- 
tions as well as modeling are needed to determine the 
luminosity distribution and decide the source. High res- 
olution imaging of submillimeter emission is a promising 
way to accomplish that. The OH megamaser in the west 
nucleus is off the peak of dust continuum and luminosity, 
and is hence unlikely to mark the position of a dominant 
heating source such as an AGN. 



We are grateful to people in and around the SMA 
project for making our challenging observations possible. 
We also thank Satoki Matsushita for the CO(2-l) data 
in Fig. [T]and the referee for comments that improved the 
paper. This research made use of the NASA/IPAC Ex- 



tragalactic Database (NED), the NASA's Astrophysics 
Data System (ADS), and the computer system in the 
Astronomical Data Analysis Center (ADAC) of the Na- 
tional Astronomical Observatory of Japan. 
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TABLE 1 
Arp 220 parameters 



parameter 


value 


ref. 


luminosity distance [Mpc] 


77.4 


1 


angular-size distance [Mpc] 


74.5 


1 


scale 


1"=361 pc 


1 


V sys (radio, LSR) [km s" 1 ] 


~5350 


2 


log(L 8 _iooo um/Lp)) 


12.2 


3 



Ref ere nces. — 1. NED; 2. IWiedner et al.l 
p002Tl : 3. ISanders et all | |20031) 



TABLE 2 

SMA OBSERVATION PARAMETERS 



parameter value 



position (J2000.0) a= 15 h 34 m 57!26 

<5=+23°30'll"4 
field of view (LSB/USB) a 35" / 34" 

center frequency (LSB/USB) 339.58 / 349.58 GHz 

bandwidth b 2.0 GHz (1732 km s" 1 ) 

spectral resolution c 3.25 MHz (2.8 km s — 1 ) 

gain calibrators J1613+342, J1635+381 

primary flux calibrator Uranus 

projected baselines for Arp 220 10.2 m - 508.1 m 

a Full width at half maximum of the primary beam in each 
sideband. b Bandwidth in each sideband. The velocity cov- 
erage for the CO(3-2) line is in the parenthesis. The ve- 
locity resolution for the CO(3-2) line is in the parenthesis. 



TABLE 3 
Log of SMA observations 



No. 
(1) 


UT date 
(2) 


configuration 
(3) 


# Ant. 
(4) 


T225 

(5) 


(Tiys.DSB) 

[K] 
(6) 


-^baseline 

[m] 
(7) 


^~obs 

[hr] 
(8) 


c*J1613+342 
°345 GHz 

(9) 


c.11635+381 
°34B GHz 

M 

(10) 


3 


2004-06-11 


1,8,9,11,14,17,18 


7 


0.07 


335 


25.8-352.9 


4.6 




(2.41) 


5 


2004-07-08 


1,5,7,8,9,11,16,23 


8 


0.08 


348 


13.5-122.9 


3.5 


1.05 


2.10 


6 


2005-05-22 


1,11,17,18,19,21 


6 


0.06 


274 


36.5-508.1 


3.4 


1.15 


0.99 


7 


2005-08-12 


1,4,5,7,8,11,23 


7 


0.08 


324 


10.2- 68.7 


2.4 


1.33 


0.71 



Note. — (1) Track n umber, (2) UT data in the yyyy-mm-dd format. (3) Array configuration. Antenna locations are 
listed fsee lHo et"aT1l2004l for a map with the numeric keys). (4) Number of participating antennas. (5) Zenith opacity at 
225 GHz measured at the Caltcch Submillimctcr Observatory next to the SMA. (6) Median system temperature toward Arp 
220. (7) Range of projected baseline length for the galaxy. (8) Total integration time on the galaxy (not necessarily with all 
the participating antennas, because of data flagging). (9) and (10) Flux densities of the quasars used for gain calibration. 
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TABLE 4 
CO and Continuum Parameters 



parameter CO (3— 2) line continuum unit 





low res. 


high res. 


low res. 


high res. 




resolution (FWHM) 


0.51x0.49 


0.38x0.28 


0.51x0.48 


0.25x0.21 


2 

arcsec 


linear resolution 


180 


120 


180 


80 


pc 


center frequency 


339.6 


339.6 


344.6 


349.6 


GHz 


velocity resolution 


30 


30 






km s -1 


bandwidth 






1.5 a 


2 


GHz 


scaling factor to Tt, 


42.4 


99.0 


42.1 


189 


K Rj/(Jy beam -1 ) 


r.m.s. noise 


38 


51 


9.8 


11.6 


mjy beam -1 




1.6 


5.0 


0.41 


2.2 


Krj 


line flux in the central 3" 


2.2 x 10 3 


1.6 x 10 3 






Jy km s 1 


continuum in the central 3" 






0.78 


0.72 


Jy 


peak integrated intensity 


1.8, 2.0 


0.8, 1.1 






10 2 Jy beam 1 km s 




7.8, 8.5 


7.8, 11 






10 3 K km s- 1 


peak intensity 


0.98 


0.46, 0.50 


0.18, 0.35 


0.11, 0.24 


Jy beam -1 




42 


46, 49 


7, 15 


20, 46 


Krj 




49 


53, 57 


14, 22 


27, 54 


K (Planck) 


CO(3-2) to (2-1) ratio b 


1.1, 1.0 










E, W line flux c 


3.0, 3.2 


2.8, 3.7 






10 2 Jy km s -1 


E, W continuum flux density* 3 






0.20, 0.38 


0.19, 0.36 


Jy 


deconvolved size d (E) 






0.31x0.25, 41 


0.27x0.14, 67 


arcsec 2 , deg 


deconvolved size d (W) 






0.28x0.22,144 


0.16x0.13,109 


arcsec 2 , deg 


deconvolved peak intensity 






43, 91 


52, 162 


K (Planck) 


spectral index q(230,345) c 






2.6±0.6, 2.7±0.5 






line width at 20% of maximum 


540, 600 








km s — 1 


kinematical major axis p. a. 




~45, ~290 






deg 


V/R 




2, 6 






km s _1 pc — 1 


R 




11, 4 






10 pc 


mean mass density inside R 




2 x 10 2 , 2 x 10 3 






M Q pc -3 


Af d y„(< R) 




11, 6 






10 8 M Q 






3, 1 






Myr 



Note. — The first value is for the east nucleus (or nuclear disk) and the second for the west when a column has 

two values separated with a comma. Flux calibration uncertainties are 15%. 

a Each sideband has 0.73 GHz. b Ratio of CO integrated brightness temperatures. The CO(2-l) data are from SS99. 
The minimum baseline and resolution are matched in both lines at 41 kA and 0'.'57, respectively. The uncertainty of 
the ratio is 18%. c Integrated within 0'.'45 from each nucleus. d The major and minor axes in FHWM and the position 
angle of the major axis from Gaussian fitting. These deconvolved sizes still include the seeing disk, whose FWHM is 
estimated to be 0'.'09 (3372j). c Spectral index, defined with S v oc u a , between 230 GHz and 345 GHz. See 



TABLE 5 
Continuum Positions in J2000.0 



parameter value 

Arp 220 E, a 15 h 34 m 57!285 
6 +23 o 30Tr!'27 
Arp 220 W, a 15 h 34 m 57!215 
S +23°30Tr:'45 
E -W offset (Act, A<5) = (0"96, -C.T8) 
E -W offset Q','98, p.a.=101° 

Note. — See Sj32]and ^5T]for astro- 
metric errors. 
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TABLE 6 

Opacity toward the nuclei 





T 860 


T 1300 


(1) 


(2) (3) 


(4) (5) 



^860/ 1300 (dust, W, 0"5)= 1.40 ±0.30 : 

2 2.3 (1.2-5.3) 1.0 (0.5-2.3) 

1 0.4 (< 3.0) 0.3 (< 2.0) 

#860/1300 ( dust . W, 0'.'23-0'.'3)= 1.57 ±0.32 : 

2 1.7 (0.8-3.5) 0.7 (0.3-1.5) 

1 ■ • • a (< 1.4) ■ • • a (< 0.9) 

#860/1300 (dust, E, 0'.'5)= 1.34 ±0.36 : 

2 2.8 (> 1.3) 1.2 (> 0.6) 
1 0.8 ( ■■■ b ) 0-6 ( ■■■ b ) 

Note. - (1) Power-law index of dust 
emissivity (oc \~^). (2) Formal solution of 
860 /J,m optical depth from Eq. ((T}. (3) 
Range of rs60 corresponding to the Ha un- 
certainty of i?860/i300- (4) and (5) Same as 
(2) and (3), respectively, but for 1300 (an 
optical depth. See £15.41 for detail. 
No solution is found. b Any value is possi- 



ble. 



TABLE 7 

Upper limit to Dust temperature Tg6o in the west 

NUCLEUS 



T 25,fg 


025 


= 0860 


025 = 0'.'25 


25 = 0"39 


T>5 


= 7860 


T 25 = 0.76T 860 


T 25 = 0.61T860 


(1) 




(2) 


(3) 


(4) 


1 


1.3 


x 10 2 


1.4 x 10 2 


1.5 x 10 2 


2 


1.8 


x 10 2 


1.8 x 10 2 


1.8 x 10 2 



Note. — Upper limits to the temperature of the 
860- /^m emitting dust (Ts6o) in the west nucleus, cal- 
culated in various possible cases. The unit is kelvin. 
See i|5.4l for detail. (1) Foreground extinction toward 
the nucleus in 24.5 /im. (2) For a less likely case in 
which the photosphere is the same for 25 ^m and 860 
fim. (3) and (4) For cases in which the 25 fim photo- 
sphere of size 825 is heated exclusively by the compact 
860 fim core of size #860 ~ 0'.'15. The t wo 25 /im sizes 
correspond to the range set by 1SN99I . The tempera- 
ture upper-limits will be lower if there is additional 
heating source between the 25 ^m photosphere and 
the 860 (im core. 
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4000 



4500 



5000 5500 
Velocity (radio, LSR) [km/s] 



6000 



6500 



Fig. 1. — CO(2-l) spectrum of Arp 220 obtained from the SMA observations of Matsushita ct al. (2008). The data cube from which the 
spectrum is obtained has a velocity resolution of 20 km s~ x and is convolved to a resolution of 5" (FWHM). Continuum is not subtracted. 
The horizontal dotted line shows the mean intensity of the outermost 40 channels of the spectrum (i.e., 20 from each edge). We adopt, on 
the basis of this spectrum, the channels outside of the velocity range of 4800—5900 km s -1 as continuum channels. The CO spectrum has a 
peak signal-to-noise ratio of 147 and sets an upper limit to high-velocity CO emission — such emission in the adopted continuum channels 
has no more than a few percent of the peak line intensity of the galaxy at this resolution. Also plotted is the (fiat) spectrum of the quasar 
J1751+096 observed during the same track and was processed in the same way as Arp 220. It shows the flatness of the line-free passband. 




349.0 349.5 350.0 350.5 

Frequency [GHz] 




4600 4800 5000 5200 5400 5600 5800 6000 
Velocity (radio, LSR) [km/s] 



Fig. 2. — (top) A spectrum of Arp 220 in the central 3" of the merger, taken from our USB data at 0"5 resolution. No continuum 
subtraction is made. The inset velocity axis is for the HCO+(4-3) line, (bottom) CO (3-2) spectrum of Arp 220 in the central 3" of the 
merger, taken from our LSB data cube of Of/5 resolution. Continuum has been subtracted. The line has a total flux of 2.2 X 10 3 Jy km s — 1 
in the velocity range of 4800-5900 km s —1 . 
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1 -1 
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Fig. 3. — SMA 0'.'5 resolution images of the CO(3-2) line and 345 GHz continuum in Arp 220. (a) CO(3-2) integrated intensity (i.e., the 
0th moment). The nth contour is at 5 X n. 1 ' 6 Jy beam -1 km s . (b) CO(3— 2) mean velocity (i.e., the 1st moment). Contours are in steps 
of 50 km s _1 . (c) CO(3— 2) velocity dispersion (i.e., the 2nd moment). Contours are in steps of 15 km s _1 . (d) CO(3-2) peak brightness 
temperature calculated using the Rayleigh- Jeans approximation. Contours are in steps of 4 K starting from 7 K. (e) 345 GHz continuum. 
The nth contour is at 25 X ra 1 15 mjy beam -1 . The first contour is at 2.5 <r and negative contours are dashed. The two plus signs are at 
the continuum peaks whose positions are in Table [3] The positional offsets are measured from the phase tracking center listed in Table |2"1 
At the bottom-left corner of each panel is the FWHM size of the synthesized beam. 



20 



SAKAMOTO et al. 



o 

CO 



CD 



O 
CD 
Q 



o 

CO 



CD 

c/> 



o 
Q 




5100 5200 5300 5400 5500 5600 



300 pc 

, , e 



-1 
R.A. offset [arcsec] 



o 

CD 

Q 




-1 - 



-2L 



0.38" x 0.28 



-1 
R.A. offset [arcsec] 

50 100 150 



300 pc 



200 



^ 



-1 



-2 



Arp 220 
0.86 mm 



Q 




Q 



0.25" x 0.21' 
1 



300 pc 





R.A. offset [arcsec] 



Fig. 4. — High resolution SMA images of the CO(3-2) line and 350 GHz continuum in Arp 220. (a) CO(3— 2) integrated intensity (i.e.. 



the 0th moment). The rath contour is at 5 



Jy beam" 1 km s" 1 . (b) CO(3-2) mean velocity (i.e., the 1st moment). Contours are in 



steps of 50 km s .(c) CO(3-2) peak brightness temperature calculated with the Rayleigh- Jeans approximation. Contours are in steps of 
5 K starting from 15 K. (d) 350 GHz continuum. Contours are in steps of 30 mjy beam -1 (=3cr). Negative contours are dashed. The two 
plus signs are at the continuum peaks whose positions are in Table [S] The positional offsets are measured from the phase tracking center 
in Table [2] At the bottom- left corner of each panel is the FHWM size of the synthesized beam. 
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Fig. 5.— CO spectra of Ar p 220 E (left) and Arp 220 W (right). The CO(3-2) spectra are from our Of/5 resolution data and the CO(2-l) 
from the 0'.'57x0''52 data of ISS99I . Continuum has been subtracted. Brightness temperatures are from the Rayleigh-Jeans approximation. 
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Fig. 6. — Positions-velocity diagrams of CO(3-2) emission across Arp 220 E (left) and Arp 220 W (right). Our 0"33 resolution data 
with continuum subtraction are used. Intensities are in the Rayleigh-Jeans brightness temperature Krj. Contours are in 1.5 <r steps with 
zero contours omitted and negative contours dashed. The vertical dotted line is at each nucleus (see Table 151 for the coordinates). The 
horizontal dotted line at 4350 km s — 1 is the approximate systemic velocity of Arp 220. The black rectangle at the bottom right corner 
shows the spatial and velocity resolutions. The spatial resolution along the cut is 0'.'36 for E and 0"28 for W. 
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Fig. 7. — CO(3-2) LTE model calculation, (a) Optical depth, (b) excess brightness temperature, and (c) equivalent width with respect 
to thermal dust emission. See text for the model assumptions. The thick contours in the plot are to guide eyes and are for tco(3— 2) = lj 

T°°^f = 30 K RJ , and CO(3-2) equivalent width of 10 3 km s _1 , respectively. The dust opacity r d is 1 at N H2 = 10 25 5 cm -2 for the 

parameters used. The panels (b) and (c) use an arbitrary Ay of 100 km s , but its choice has little effect on the estimate of Nn 2 /AV 
in Arp 220 (see Appendix [Bj) . 
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Fig. 8. — Compact cm- wave sources and sub-mm continuum emission registered with it. Red plus signs are compact (< 1 pc) continuum 
sources thought to be radio supernovae or young (30-50 yr) supernova remnants (Lonsdale ct al. 2006; Parra ct al. 2007). The red circle 
in each nucleus indic ates the 'center' (= median position) of the compact sources. The blue circles and ellipses are OH megamasers 
(Lonsdale ct al. 1998). The green circle is the peak of diffuse cm emission of the west nucleus, as measured with respect to the compact 
sources bv lRovilos et al.l $2003 ). The contour map is from our 860 /an continuum data, and is shifted so that the west peak coincides with 
the 'center' of the supernovae distribution in the nucleus. The black squ ares show the positions of the sub-mm peaks before the shift. The 
coordinates of the black squares and red circles are in Table \E\ and £|5.1I respectively. 




delta-R.A. [arcsec] delta-R.A. [arcsec] 

Fig. 9. — Spatially smoothed distri bution of the s up ernovae and supern ova remnants in Arp 220. Each plus sign is a compact (< 1 pc) 
cm continuum source listed by Lonsdale et al. (2006) or Parra ct al. (2007) as a radio supernova or a young (30—50 yr) supernova remnant. 
The contour map is made by convolving the point sources with a Gaussian kernel whose FWHM is shown at the bottom-left corner. All 
sources are given the same weight to show smoothed spatial distribution of source density and of massive star formation rate. Contours 
are in 10% steps starting from 10% of the maximum. 
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APPENDIX 
1. DISK-GAUSSIAN DECONVOLUTION 



One sometimes needs to estimate the size of a uniform-brightness source from its marginally resolved map. We 
here note that the diameter d of such a disk source is about 1.6 times the FWHM one would obtain from the usual 
Gaussian-source deconvolution, in which the well-known quadrature relation holds 

Cap = ^int + OpSFi (Al) 

where map , #mt, and 6psf are the FWHM size of the source in the map, the intrinsic FWHM size of the Gaussian 
source, and the FWHM of the Gaussian point spread function, respectively. 
A uniform-brightness circular source of the diameter d and total flux 1 has the visibility 

where p is the baseline length and J is the Bessel function of the first kind. A Gaussian source of FWHM 0; nt and the 
same total flux has the visibility 



T 2 e 2 c 2 

4 In 2 

The half-power uv radius of the visibility functions is at ndp ~ 2.2151 for the disk and Tv9- m tP = 2 In 2 for the Gaussian. 
The two functions agree at their half power points if 

2.2151 

dw — r- — int « 1.60 6 int . (A4) 
2 In 2 

Fig. [TU]shows that the two visibility functions resemble each other out to longer baselines when the diameter is slightly 
smaller than this. (The figure is for d = 1.57 6*i nt .) 
Taking the convolution of the disk with the PSF more directly, the convolution relation in the uv domain is 

WI^^W). (A5) 

P I 41n2 J ndp P \ 41n2 J V ' 

Requiring this to be exact at the half-power point (i.e., at p that satisfies 7r# map /9 = 2 In 2) one obtains 

C p = ^ 2 + ^sf (A6) 



for d/9 map <C 1. Comparison of this with I All gives a similar coefficient as in lA4[ 

d« 1.70 flhrt. (A7) 
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Fig. 10. — Visibility amplitudes of a Gaussian of 6*fwhm = 0'.'15 (dashed line) and a disk with a diameter d = 1.570pwHM (solid line). 
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2. EQUIVALENT WIDTH OF MOLECULAR LINE WITH RESPECT TO DUST CONTINUUM 

The equivalent width of a molecular line to the dust continuum emission for an isothermal mixture of molecular gas 
and dust has the following relation to the ISM properties. Let r m and r d be the optical depths of the molecular line 
and the dust continuum adjacent to the line's frequency is, respectively. The total optical depth at the line frequency is 
T"m + 7d- The emission intensities on and off the line frequency are, for LTE excitation at the gas (and dust) temperature 
T, 



I» ne+cont -=B v (T) l- e -^m+T d ; + B ^ TcMB y-(T m +r d) ( B1) 

I^ =B V {T) [1 - e- T d] + B,(T C MB)e- Td , (B2) 

where B v is the Planck function and Tcmb is the cosmic background temperature. The excess brightness temperature 
of the dust continuum is the excess of intensity, expressed in the Rayleigh- Jeans temperature, over an off-source 
position. 

^bTxccss = ^ [C nt - M - C nt ' (0)] = (1 - e- Td ) [J(v, T) - J{v, T CMB )] , (B3) 

where J{y, T) = 2j£(e hv ' kT — l) -1 is the effective radiation temperature and c,fc, and h are the speed of light, Boltzmann, 
and Planck constants, respectively. The excess brightness temperature of the line after continuum subtraction is, 

^Scess = ^2 { [4 inC+C ° nt -(^ + T d ) - /i inC+C ° nt -(0)] ^ [C^Td) - C^O)] } 

- e- Td (1 - e- T -) [J{v, T) - J{v, T C mb)] • (B4) 

Hence, the equivalent width for a line with a velocity width of AV, is 

T iino A y 1 _ Tm 
b,cxcc BB = i_e 

e q /yncont. „T d 1 V ' 

b, excess 

The beam filling factor of the ISM, which should be multiplied to Tbcxccss an d ^b excess to compare them with obser- 
vations, cancels out in W eq . 
The line optical depth of a CO-like molecule in LTE is written in the following form 

r m ~ a^T-" (B6) 

at high temperatures (T 3> hv/k), and the dust optical depth can be written as 

T d = WV H2 , (B7) 

where A/h 2 is the H2 column density, a is a constant that reflects the properties of the molecule in question and its 
abundance with respect to H2, and b is another constant that contains gas-to-dust mass ratio and the dust mass opacity 
coefficient. The equivalent width can be uniquely calculated on the Nn 2 /AV - T plane when the dust emission is 
optically thin (i.e. e Td -1» Td). Moreover, if r d <C 1 then W cq depends on only one of the two parameters except 
when the molecular line has a moderate opacity. 

w _ / (bNnJAV)- 1 (for r d < 1 and r m > 1) , , 

^ cq ~ \ ab^T- 2 (for r d < 1 and r m « 1) 

This is why the iso-We q contours are mostly parallel to either the Nb_ 2 /AV or T axis in Fig. [7] (c) and why W cq 
can be used to estimate one of the two parameters. The N-a 2 /AV estimated using the approximate formula (|B8p 
overestimates the correct one from (|B5|) by a factor of / ln(l + ^-), or by less than 14% for AV^ < 300 km s _1 and 

W cq = 10 3 km s _1 . Finally, the factors that can bias the above evaluation include non-LTE excitation, disagreement 
of gas and dust temperatures, non-uniform ISM density, column density, and temperature, and variations of molecular 
abundance and gas-to-dust ratio. 

3. DUST CONTINUUM INTENSITY AND SN SURFACE DENSITY 

Surface density of supernovae and supernova remnants (called SN in short), Ssn, must be approximately proportional 
to the intensity of dust continuum, 7 d , for a dusty starburst region if some plausible assumptions hold. First, if the 
number of SN (-ZVsn) i n a region is proportional to the star formation rate Rsf in the region then we have 

iVsN°C.RsF (CI) 

S sn cxS S f- (C2) 

The second formula expresses the first one in terms of surface density. Secondly, assuming that the region is dominantly 
heated by star formation (of luminosity Lsf) and cooled by dust emission, we have 

XsfocT 4 / 2 (C3) 
.-. SsfkT 4 , (C4) 
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where / and T are the size and temperature of the region, respectively. Thirdly, assuming that star formation surface 
density is related to the gas surface density £ g , we employ the Schmidt-Kennicutt law 

SsF^Eg™, (C5) 

where n is widely thought to be in the range of 1-2 (|Kennicuttlll998| ). Finally, if the dust emission that we observe in 
(sub)millimeter is of moderate optical depth or optically thin (i.e., < 1 for dust opacity) then we have 

Id oc TdT oc E g T. (C6) 

It follows from these relations that 
and 

/ d cxS SN " +i . (C8) 
The exponent — + i is in the range of 0.75-1.25 for n =1-2 and is 0.96 for the canonical n = 1.4. 



T«E SN 1/4 (C7) 



